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We generalize the effective field theory of single clock inflation to include dissipative effects. 
Working in unitary gauge we couple a set of composite operators, 0^,^..., in the effective ac- 
tion which is constrained solely by invariance under time-dependent spatial diffeomorphisms. 
We restrict ourselves to situations where the degrees of freedom responsible for dissipation do 
not contribute to the density perturbations at late time. The dynamics of the perturbations 
is then modified by the appearance of 'friction' and noise terms, and assuming certain local- 
ity properties for the Green's functions of these composite operators, we show that there is a 
regime characterized by a large friction term 7 ^ _ff in which the (^-correlators are dominated 
by the noise and the power spectrum can be significantly enhanced. We also compute the 
three point function (CCC) ^ wide class of models and discuss under which circumstances 
large friction leads to an increased level of non-Gaussianities. In particular, under our as- 
sumptions, we show that strong dissipation together with the required non-linear realization 
of the symmetries implies |/nl| ^ -^rj^ ^1. As a paradigmatic example we work out a 
variation of the 'trapped infiation' scenario with local response functions and perform the 
matching with our effective theory. A detection of the generic type of signatures that result 
from incorporating dissipative effects during infiation, as we describe here, would teach us 
about the dynamics of the early universe and also extend the parameter space of inflationary 
models. 
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INTRODUCTION & MAIN RESULTS 



The Effective Field Theory (EFT) paradigm is one of the cornerstones of theoretical physics, 
from the standard model to condensed matter systems [H, Q]- EFT ideas have recently gathered 
thrust also in the realm of gravitational physics. For example, EFT techniques have been introduced 
in to solve for the dynamics of coalescing binary systems to great accuracy [fi-ll]; and an 



EFT setup has been proposed for the study of cosmological perturbations in |l2l |. 

The EFT of inflation for the case of single field (one clock) models was developed in 
The starting point is an action in unitary gauge (where all the fluctuating degrees of freedom 
are encoded in the metric) which is required solely to be invariant under time dependent spatial 
diffeomorphisms. The advantage of this approach is that it enables us to parameterize all possible 
signatures of inflation in terms of a set of coefficients for ('higher-dimensional') operators in a 
Lagrangian built with the low energy (large distance) degrees of freedom, and constrained only by 
the symmetries of the theory. Within the EFT it is possible to describe the fluctuations around 
an approximate de Sitter background without any assumption about the fundamental degree of 
freedom that is driving inflation^ (which may as well be strongly coupled). To that end it is useful 
to restore time diffeomorphisms (broken by the existence of a preferred time slicing) by means 
of the Stiickelberg field tt, which is the Goldstone boson that realizes time reparameterizations 
non-linearly. 

There are two important consequences of introducing the vr field which turn out to be extremely 
helpful. First of all, one notices that at sufficiently large energies {E » \^H) the Goldstone boson 
captures all the information about the physical scalar mode (C — —Hit), namely the 'longitudinal 
mode' in unitary gauge. ^ And secondly, there is a limit (decoupling limit, i.e. e — )■ 0) in which 
we may ignore all the effects induced by gravitational interactions.^ These two observations allow 
us to concentrate on a theory of Goldstone bosons, whose interactions are dictated by symmetry, 
which greatly simplifies the computations. 



One may wonder about the underlying theory that produces the background. However, once the inflationary 

paradigm is accepted, it is ultimately the theory of fluctuations that is directly tested by observations. 
^ This is similar to what occurs in gauge theories, for example in longitudinal WW (or Z) scattering, whose amplitude 

can be obtained in terms of processes involving the associated Goldstone bosons ('eaten' by the W's and Z's) at 

high enough energies, E ^ mw ('equivalence theorem') [3, [2^. 
^ This is also equivalent to sending the gauge coupling g to zero while keeping the symmetry breaking scale v fixed, 

or in our case taking Mp — >■ 00 and keeping Mp\H\ finite. 
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An EFT for multi-field inflation was recently introduced in [2g], where new (light) degrees of 
freedom were included. Multi-field inflation can reproduce the signals from single field models, but 
can also give rise to new ones which (provided certain shapes are detected) may allow us not only to 
separate between the two, but also distinguish amongst different realization in multi-field scenarios. 
In a the authors concentrated on the case in which additional degrees of freedom (ADOF) affect 
directly the overall curvature perturbation (or isocurvature perturbations). This can happen for 
example if these extra fields modify the reheating time and consequently the duration of inflation, 
or if they affect the composition of the plasma in the reheating phase. Because of this, the fields 



considered in 26|] were light scalars so that they acquired scale invariant perturbations. 



In this paper we also consider situations in which ADOF other than the Goldstone boson are 



excited. However, contrary to the models in 26|], we concentrate on cases where this extra sector 



does not directly affect the duration of inflation, or the composition of the plasma, but it alters 
the dynamics of inflation by directly coupling to the clock around or before the time the modes we 
observe cross the horizon. This includes, for instance, the 'trapped inflation' scenario [27*] where the 
production of particles modifies the evolution of the inflaton <j), while producing negligible direct 
contributions to (late time) density fluctuations due to dilution. 

Since in general these new particles will contribute to the stress energy tensor of the background, 
the true Goldstone boson includes not only the physical clock, which we assume uniquely controls 
the physics of inflation, but also a component that depends on the fluctuations of these ADOF. 
To isolate the relevant component whose perturbations control the observed density fluctuations, 
here we will not work with this Goldstone boson (namely the field whose quadratic Lagrangian 
is uniquely fixed by the background), and thus reserve vr for the fluctuations of the clock that 
determines the end of inflation (for example, the inflaton cf) in the model of [23] )• This choice 
will slightly modify the construction of the EFT, and in particular the choice of 'unitary gauge' 
and overall normalization of our vr. However, it maintains the relationship between vr and i.e. 
~ —Htt (at linear order), which we find more convenient. 

Given that the ADOF do not explicitly contribute to observable quantities, this suggests we 
may integrate them out and obtain an effective action in terms only of vr. However, this procedure 
is not straightforward, mainly for two reasons. First of all, even though we do not observe the 
fluctuations of the ADOF as 'external states', they are produced during inflation and in general 
are not in the vacuum; and secondly, these fluctuations may not be gapped, or in other words, we 
are allowing for very soft (essentially gapless) collective excitations. This means that the effective 
action cannot be described with only one degree of freedom, and in turn this will be linked to 
dissipation whose effective description is one of the goals of the present paper. Notice that this 



does not mean that the ADOF are necessarily 'light'. For example, as shown in [27|, heavy 
particles (compared to the Hubble scale) can be ultimately produced by the time dependence in 
the Hamiltonian induced by the physical clock (j).'^ However, their influence in the dynamics of the 
perturbations of the clock, which includes dissipation and noise, remains active even at low(er) 
frequencies, i.e. a; ~ if. 



One can study dissipative effects using the 'in-in' closed-path-time formalism 29(], however, 
here we will resort to a different setup. Similarly to the EFT for dissipation introduced in 13, 31] 
(to deal with gravitational wave absorption in binary black hole systems), we will incorporate 
dissipative effects in the EFT of inflation by coupling the metric in unitary gauge to a set of (scalar, 
vector and tensor) composite operators, O^j^..., constrained solely by the symmetries of the EFT. 



In the case of trapped inflation particles are produced when the adiabaticity condition is violated [27I. [2^. 
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All the information about the dissipative sector is thus encoded in a set of correlation functions 
which can be matched against observation, e.g. the power spectrum, non-Gaussianities, etc. (A 
similar formalism [s^] can be used to describe dissipative effects in the EFT for hydrodynamics 
developed in [s^.) This approach is clearly very ample, and for that reason it is also difficult to 
treat in complete generality. However, it is possible to identify a physical regime in parameter 
space where many simplifications occur. We will consider physical situations where the time scale 
for dissipation and fluctuation induced by the ADOF is much smaller than a Hubble time, with 
negligible memory effects. Under this condition fluctuation and dissipative effects become local, 
which allows us to study many possible realizations in complete generality. We spell out our 
assumptions in detail throughout the paper. 

Two separate type of contributions enter in the dynamics of vr. There is the noise induced by 
the ADOF, and there is also the back reaction due to the mutual interaction, i.e. 'friction', namely 
a 771" term. (This is due to the physical fact that during inflation energy is being damped into the 
ADOF.) Here we concentrate in the strong dissipative regime, where 7 is taken to be much larger 
than H. In such circumstances, we show the memory on the initial conditions washes out and 
the power spectrum and non-Gaussianities are dominated by the noise induced by the ADOF. 
Moreover, the former can be significantly enhanced with respect to the quantum fluctuations in 
the Bunch-Davies state. The reason is twofold: first of all the size of the fluctuations for the 
normalized tt field (tTc = -v/iVcVr) is larger than Hubble, this is because these are not produced by 
the vacuum (also they freeze out at a larger value of k/a^ i.e. Csk/a^ = Csk^ — ^ fo'^ 

7 ^ H); and secondly, the normalization scale Nc can be smaller than what is required by Einstein 
equations in the absence of any other contribution to the stress energy tensor, i.e. Nc < 2Mp\H\/c^. 



In the EFT of inflation [13[] the non-linear realization of time diffeomorphisms allows us to 
relate different observables (such as the two and three point ^-correlation functions). For instance, 
considering the case of scalar couplings and taking the noise to be Gaussian, we show that terms 
of the form S^^oo = / ^J—gOg^^ (in unitary gauge) can increase the level of non-Gaussianities by 
a factor oi ^/H with respect to the result for single field inflation without ADOF [l^, yielding 



I/nl I — l/{HCs). In the spirit of the EFT for multifleld inflation of [26| we discuss a class of 
models which flt into this category. On the other hand, operators of the form Sfo = / yj—gf{t)0 
produce two types of non-linearities. Either sourced by direct non-linear couplings to O, or induced 
by contributions beyond linear response theory. We show there is a class of models where the 
computation of dissipation and non-Gaussianities are linked. Assuming as we do the existence of 
a preferred clock that drives inflation, we show that the linear and non-linear response are indeed 
related in such a way to produce non-Gaussianities of order |/nl| — l/ic^H). 

The basic idea is perhaps better illustrated in the standard inflationary scenario of a slowly 
rolling scalar field (j). As we know from our classical mechanics intuition, to induce dissipation we 
need factors of (p in the equation of motion (EOM). On the other hand, (non-linear ly realized) 
general covariance requires dt(f) — )• n^^df^tp where is the normal vector orthogonal to the equal 
time surfaces given by n'^ = g^^ dyC^j \J—{d(^ . At the level of the perturbations of the clock, 
namely — )• -|- this induces a dissipative term but also non-linear interactions, and in 
particular a term ^{didcj))'^ /(j) (properly normalized). The latter gives non-Gaussianities of order 

lidMf f A ^ If I 7 /-.s 
' /nlC I/nlI ~ (ij 



c 



where we used C, ~ —H6(j)/(j), and allowed for Cg < 1. We will make this argument more precise 
and at the same time generic for all models of single clock inflation. We conclude that in either 
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case, 50^00 or Sf(i-^Q, the non-linear interactions are significantly enhanced in the strong dissipative 
regime with 7 3> -ff, as one would have naively expected. 

Finally, vector couplings such as j \/—g0^g'^^, may induce a large friction term but without 
the addition of non-linear couplings between vr and C^. However, once again depending on the 
model, the non-linear response will produce large non-linearities as above. 

This paper is organized as follows: In the remaining of sec. [T] we discuss the basic ideas and 
results, putting emphasis on the overall picture rather than the technicalities of the calculations. 
Then in sees. [2] -[7] we explicitly construct the EFT to include dissipation in inflation and provide 
detailed support for our claims. In sec. [8] we perform the matching for a key example: (a local 
version of) trapped inflation. The idea of including dissipative effects during inflation is also a 



key element of the warm inflation paradigm 34l436f|. which we will briefly comment upon towards 



the end. We relegate other examples and more technical points to appendices. Everywhere we set 
c = h = ks = I and adopt the mostly plus sign convention. 

1.1. Preliminaries 

Let us imagine adding a friction term to the (one-dimensional) harmonic oscillator. For conve- 
nience of notation let us denote the displacement from equilibrium as iT{t). The EOM reads 

vi" + 7^ + ^0'^ = J, (2) 

where for future purposes we added a stochastic force with (J) = 0. As it is well known, this 
seemingly innocuous equation does not derive from a Lagrangian of the form £(vr, tt).^ The reason 
is simple, energy is not conserved. In fact. 

The expression in Eq. ([2]) is local in time, however, in general the effective EOM for vr takes 
the so called 'Langevin' form, which is non-local, i.e. 

jf + UjlTr+ [ dt'^it - t')TT{t') = J{t). (4) 



Clearly, our treatment would greatly simplify if we were allowed to perform a local approxima- 
tion. This in turn amounts to making an assumption about the ADOF leading to ^{t) in the above 
equation. In particular, to get the EOM in the form of Eq. ([2|) we need 

Im7(a;) ~ 70;, (5) 

with 7 a constant. The relationship in Eq. ^ is sometimes referred in the literature as Ohmic 



behavior [371]. We illustrate some examples in appendix Rl 



In practice we do not expect Eq. ([5]) to hold up to arbitrarily high frequencies, and a more 
realistic (though phenomenological) approach is provided by Drude's model, with Eq. ([5]) replaced 



^ One can, nonetheless, construct models where a dissipative equation results from an effective description where 
the ADOF responsible for dissipation are 'integrated out' in a 'in-in' formalism with twice as many degrees of 



r esp 

freedom [29l. IsTI. l38l|. This is not the route we follow in this paper (see next sub-section). 
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by [37| 

Im7z)(a;) = 7W (l + (^"^ /uj^,) ^ , (6) 

where ujd serves as a cutoff. In this case there is a memory time on the scale r^) ~ , and an 
exponential damping 

71) (t) ~ e(t)7D . (7) 

Instead of relying on assumptions about the physics of the ADOF, sometimes it is possible to 
connect the properties of the noise and Green's functions by means of some well known theorems; 
for example if we assume the noise satisfies the following conditions: 

{Jit)) = (8) 
{J{t)J{t')) ^ vj5{t-t'), (9) 

with vj a constant. This is the case, for instance, if the "environment" (i.e. the dissipative ADOF) 
is placed at afsufficiently large) temperature T. Then, using the Fluctuation-Dissipation (FD) 
theorem 2^, 37 1 one can show 



Im7(w) ~ (10) 

(assuming equilibration occurs sufficiently fast after the perturbations are turned on), or in other 
words vj ~ ^T. In this scenario the EOM becomes local and memory effects are washed away. 

To keep the treatment as simple as possible, we will study situations where the local approxi- 
mation applies. Later on we will discuss some specific examples. However, at the level of the EFT 
we refrain from adopting any model for the underlying dynamics of the ADOF. We introduce the 
basic idea of our approach next. 



1.2. The story of O 

We consider now the generic situation where we have a theory for vr that describe small (long 
wavelength) perturbations of a dynamical system with Lagrangian C^^. Following ^,i31] to include 
dissipation we couple vr to a composite operator O such that 

Sn,t = - J d^xO{x)Tr{x). (11) 

For the cases where there is a shift symmetry, vr — )• vr + c, the interaction takes the form 

^int = J d^xd{x)7r{x), (12) 

which can be described as in Eq. (Ilip by replacing O — t- C 

The addition of S'int allows us to study the response of our system to the interaction with 
a dissipative sector (represented by O) in complete generality. The virtue of this approach lies 
in that we do not need to assume any specific representation for the dynamics of O (which 
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could in principle represent a strongly coupled sector), and we just need to make sure inter 



actions such as Eq. ([IT]) respect the symmetries of the long distance physics described by n [30|,l3l|- 



We define now an operator 60 = O — O, with O the background expectation value of O. We 
then split 60 into two pieces, schematically, 

50 = 50s{x) + 60r{x), (13) 

with 60s{x) representing the stochastic part of 0{x) in the absence of vr, whereas 60ji{x) corre- 
sponds to the change of the expectation value of O that results as a response through the interaction 
to a TT fluctuation. We start by computing 50ii{x) within linear response theory, with vr playing 
the role of the external 'force' that disturbs the dynamics of the degrees of freedom associated with 
the operator O. (For an introductory account of response theory see for instance f39|.) Namely, 

50Rix) = - j d^yGZt{x,y)7T{y), (14) 

where 

Gg,{x,y) = i{[60{x),60{ymt,-ty). (15) 

If we denote as Dj^tt = the linearized EOM that derives from C^^ in the absence of ADOF, 
then the addition of Eq. (ITT]) leads to 

D^vr = -C + 0(7r2). (16) 

After we solve for O we get 

D^TT - J dScZtix - y)vr(y) = -60s + 0(^'), (17) 



or in Fourier space 



D^{ci,uj) - G?,t(q,c^)) vr,(^) = -<50s(q,^) + • • ■ ■ (18) 



In the above expression /^^(q, w) is an analytic function, but presumably this may not always 
be the case for G^^(q, tj). However, we will assume we can use a local approximation for the 
dynamics induced by the Green's function for a vast range of frequencies up to corrections of order 
{q/ Mq , oj /T o) , where Mo,Tq are the typical scales at which the non- locality starts to be non- 
negligible. We will make this more precise when we study the inflationary case. At any rate, it is 
clear that in order to recover an equation as in ([2]) we need^ 

lmGg,iuj, q) ^ lmGg,iu;, 0) ^ juj, (19) 

for cj > 0, with J = —50s- While an expansion in derivatives is justified by the fact that we 
can work at a; ^ Fq, the order at which the Taylor expansion starts is an assumption about the 
UV physics. (Notice that this is the lowest analytic order since the imaginary part of the Green's 
function has to be odd in w.) 

Contrary to the imaginary part, the real part of the Green's function is even in uj, and 



^ As one would expect dissipation is associated with the imaginary part of the Green's functions in Fourier space, 
a.k.a. the optical theorem, see appendix [B] 
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therefore need not vanish as w — t- 07 This would lead to a mass for vr (and consequently for 
C), hence to an evolution for outside the horizon. This is not forbidden by any principle, but 
if this was the case it would lead to an effect on curvature perturbations at late times from 
the ADOF, which ought to be negligible by our main assumption. Therefore, in this paper 
we consider situations where the ADOF are sensitive only to derivatives of vr rather than the 
value of TT itself. This requires an effective shift symmetry at the level of the response. (This 
is the case for example in the model of |27|] where the production of ADOF is connected with 
a breaking of adiabaticity [2^, due to the time dependent background, and fluctuations in 
the ADOF are related to derivatives of vr.) If for instance inflation is driven by a scalar field 
(f), in principle the response of the O operators might be very non-linearly related to (p, and 
the Green's function for vr may as well depend non-linearly on the coupling constants in the 
Lagrangian that includes both (p and O. In fact, the shift symmetry may not be even present 
at this stage. For this reason we will employ the term 'emergent shift symmetry' to refer to 
cases when the real part of the Green's function vanishes as w — )• 0. We give an example of 



this phenomenon in sec. lS.ll where we discuss a simplified version of the trapped inflation model 27l |. 



In principle if we had a specific UV model of the ADOF in mind, we could compute the exact 
Green's function and compare with the one derived from the EFT. This is typical in EFTs and it 
is often referred as matching [Ij] . We will provide a realization of this matching procedure when we 
study some specific models later on in sec. [HI (See also appendices [C] and [Hi) Nevertheless, even 
though we do not have an explicit description for O (which might be rather involved), the key 
point of the EFT approach is that we can still study the dynamics of vr in terms of the Green's 
function of the type of Eq. (fT5|) under the approximation of Eq. (fT9l) . where 7 is kept as a free 
parameter. 

Notice that the condition for a 771" local dissipative dynamics implies 

ImGZd^) ~ (20) 

for the operator O in Eq. (jl2p . This behavior is not allowed near a; ~ by some basic analytic 
properties of the Green's functions. One can nonetheless imagine situations where Eq. pop holds 
at intermediate frequencies, i.e. fio ^ w ^ Tq (with related to the response functions of 

the Cs), while at very low frequencies the Green's function is analytic. Since we are assuming the 
relevant energy scales are smaller than the typical ones in the ADOF sector, that we take to be 
of order Fq, in order for such behavior to occur it requires the Green's function to have a mass 
scale anomalously low compared to Tq. This is (most probably) a sign of tuning in the effective 
theory.^ In any case, it is reasonable to assume this tuning would affect only one parameter (see 
appendix [H] for more details). 

An example the reader may be familiar with is the so called Abraham-Lorentz-Dirac (ADL) 
force, which arises from the (velocity dependent) interaction A-fr after we integrate out the 



^ Moreover, the real part of the Green's function G^t(^iCl) ma^y a^lso contribute to the speed of sound via a term 
quadratic in q. 

* We show in appendix iHl how, upon tuning, one may obtain the scaling of Eq. (|20p . in particular 

3 

uj'lmGZtH ^ 7 2 + O{uj/ro) ~ + O{^lo/u;,0J/^o) (21) 

for fio ^ u) ^ To, such as it is required for a 77r term stemming from the coupling O-k. 
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electromagnetic field (here A plays the role of O). The EOM turns out to be local in time [4C 
provided we choose boundary conditions where all radiation is outgoing (that is we do not include 
'mirrors').^ However, in this case the dissipative term depends on the third time derivative of the 
position, i.e. 'Tr'/Ae with Ag — rrie/e^ a cutoff scale related to the unitarity bound of the theory. 
(This follows from an expression similar to Eq. (I19p applied to the vector potential.) Even though 
in this paper we mostly concentrate on dissipative effects represented by Eq. ^ we will also 
comment on higher derivative couplings in sec. 17.51 

As we mentioned earlier another setup where a local approximation appears naturally is to 
consider a white noise for the ^O^'s, as it would be the case in thermal equilibrium at large 
temperatures. In such scenario 

(50s(k,i)(5C»5(q,t')) = {2^fyo5{t-t')5^^\]^ + ci), (22) 
and using the FD theorem we get (see Eq. (fTOl) ) 



uj ' T 

as required. 



ImGgt(w) _ vo , . 

= 7 = 777. (23) 



1.3. The two-point function 

One of the most important observables we are interested in this paper is the two-point function 
of the vr field at horizon crossing, which is related to the two-point function for the curvature 



perturbation ((" ~ —Hn 13|]), that is conserved outside the horizon. In the standard scenario of 



an expanding universe the linearized EOM for the vr field is equivalent to Eq. ([2]) with 7 — )■ 3-ff 
and ujQ = CskpYi = that is (notice that now we have a time dependent ojq) 

^2 

% + 3i?^fc + c2— vTfc = 0. (24) 

The reason the mode freezes out is due to the fact that the term proportional to ojq goes to zero 
as t — 7- -l-oo, for a fixed (co-moving) A;, and a constant value for vr solves the equation. The time at 
which this happens is determined by the condition k^, = k/a{t*) ~ H{t*)/cs, or lo^ ~ H^,. (The :*r 
denotes a quantity at freeze out.) If we impose the Bunch-Davis state as initial condition, the well 
known result is |13(]: 

with e = —H/H"^. Since these are the quantum zero-point energy fluctuations, this expression 
follows straightforwardly from 



Non-local dynamics appears in the so called 'memory effect' in gravitational wave radiation off coalescent binary 
inspirals i41]. The latter is entirely due to the non-linear interactions of the gravitational field which are not 
present in electromagnetism. 
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Naively one would expect that the addition of a friction term to Eq. (|24p . of the form 777^, will 
modify the crossing condition to k/a{t*) ~ j{t*)/cs, leading to ~ 7. Hence, from Eq. (p6|) . 
it appears as if it would produce a larger two-point function. However, this is incorrect, and the 
contribution from the homogenous equation turns out to be negligible. We will show this in detail 
later on (see sec. 15. 2p . but the basic idea is rather simple as we argue next.^'' 



1.3.1. Homogeneous solution 

Let us take Eq. ()24p but assume there is an extra dissipative term 771"^, with 'j ^ H. To gain 
some intuition we will solve the equation adiabatically starting with constant values of ojq = Csfcph- 
(For reasons that will be clear in sec. 15.21 we also parameterize time running from [— |to|,0].) This 
is a good approximation as long as wq/wq <C 1, which holds provided Cskph > H. (In fact, as we 
will see, c^A;* ~ Vt^ ^ H .) It is easy to see there are two independent solution, namely 



f^it) = exp 




(27) 



We will consider two regimes. First, at some early time (|to| ^ I/t) we assume we are in the 
solution with wq = c^feph ^ 7, so that we match it with the usual oscillatory behavior normalized 
to the Bunch-Davies vacuum 

/4 = -^e±^"«*o. (28) 



This is justified by realizing that by going sufficiently back in time, the mode begins to oscillate 
fast enough to decouple from the ADOF. We expect this to happen for lo <; Tq. This fixes the 
overall coefficient to 

-7|tol 



^± - ^7^ + 0(7M)- (29) 



As time progresses, we enter our second regime, where ujq decreased to the point the mode freezes 
out, Wo — )• Wg = Cgki,. One can then show that after matching both regimes the solution that 
dominates scales like^^ 

_2lM , *,2, 
p 2 ~l"nJ * 

f---^e^^. (30) 



As a consequence the homogenous solution acquires a damping factor e 2 ^ 1. A detailed 
analysis shows that indeed it acquires this type of exponential suppression (see section 15. 2p . 

From here we conclude that the contribution to the two-point function from the homogenous 
part becomes exponentially small as t — )• 0, which opens the door for the source noise to dominate. 



Notice that the term in 771- does not have the same role as the standard SZ/tt one in an expanding universe, since 
for the latter the frequency of all modes redshift at the same rate H , while that is not case with 7. 
The other solution decays faster with time. 
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1.3.2. Noise 

In what follows we present a basic physical argument for the computation of (vrvr) due to the 
noise dOs- (For ease of notation here we return to the more traditional range for time, t £ [0, oo].) 
The detailed analysis will be presented in section [6l 

In the EOM of ([2]) we now have to deal with an extra term, namely 

Tfk + l^k + ^iT^k = -N-^60s, (31) 

where Nc is a normalization factor. Since we take 7 ^ -ff, we work in the limit where 

ujo « 7- (32) 

Under this condition we can in principle find the exact Green's function, however it is easier to 
look at the simplified version that holds in the overdamped limit 



(33) 



where we drop the factor -^^f- ^ 1~d^- The solution reads 

G^Jt-t') = -e-'^^'''''^e{t-t'). (34) 
r 7 

(The k dependence is implicit in ujq.) We see that the response induced by the Green's function 
is approximately constant for sources concentrated on late times, while it becomes exponentially 
damped for very early sources. This allows us to define an 'equilibration time' as the scale con- 
trolling the exponential suppression: 

~ ^. (35) 



The solution for vr now reads 



poo 

TTkit) = -N-^ / dt'G'^it - t')60s(k, t'). (36) 
Jo 

(This is all what is left at late times t ^ Teq since the homogeneous solution dies away.) Assuming 
a white noise spectrum, 

{6Os0^,t')6Os{q,t)) ~ {27Tfuod{t-t')6^^^{ci + k), (37) 

the two-point function turns into 

/"OO /"OO 

{nk{t)TTg{t)) = N-^ / dt"dt'G';it-t')G^t-t"){S0s{Kt')50s{<l,t")) (38) 
Jo Jo 

~ (2^)3<5(3)(k + q)i/oiV-2 dt' (G';{t - t'i 
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Performing the integral we obtain 

(T.(*K(*)) ^^fl-.-"^') + ^ ^^f"(k + q). (39) 

which tends to a constant as t — t- +oo as we expected. 

In thermal equilibrium, when the FD theorem applies, we can relate the amplitude of the noise 
vo to the damping scale of the Green's function and the temperature T. If that was the case we 
would then have 

vo = N,^T (vr^)--^, (40) 

or equivalently 

N,ujI{tt^) ~ T. (41) 

This expression is suggestive because it reminds us of the equipartition of energy in thermal 
equilibrium. If we interpret NcUJq as a spring constant kg and tt as an harmonic oscillator, then 
ksi'Tr'^) ~ T. Indeed the factor of Nc is the canonical normalization for the field vr, and in a 
sense it represents the 'mass' of the harmonic field vr (which does not have a mass in the strict sense). 

The above expressions allow us to understand the properties of the Green's functions in the 
expanding universe, which is the case of interest here. The equilibration time represents the time 
it takes for the interactions to cancel out the effect of an initial fiuctuation. This effect is due to 
dissipation. Indeed, in the absence of dissipation the effect of the initial conditions never disappear. 
The expression in Eq. (I39p . which was obtained in a Minkowski background, is also valid in the 
limit in which we can neglect the time scale of variation of loq, given by H~^, with respect to the 
time scale of the Green's function, i.e. T~g^. When this condition is violated we cannot trust the 
solution any longer. However, we can still look at the EOM for vr and realize that since it contains 
only derivatives, if the noise is sufficiently concentrated at short distances, the correlation function 
becomes a constant. Using ujq = Cgk/a this happens when 

H Wo ~ =^ ~ a/ (42) 

Notice that at freezing the physical momentum is much larger than H for 7 S> i?. Then from 
Eq. (j39p and matching the solution deep inside the horizon and at horizon crossing we obtain 
(after re-inserting the factors of a, the scale factor) 

where we used (wq)^ = (csk/ai,)"^ and 1/a^ = y/^yH^/{c*k). For instance if we use the relation in 
Eq. (gD]) we get 

(vr,vr,)(t.):.(2vr)3^^5(3)(k + q). (44) 

These are indeed the results we find in the full computation (see Eqs. (jl32[ 1134]) ). Notice that 
depending on the value of z^o, 7 and Nc (and/or T), the two point function can be significantly 
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enhanced with respect to the standard result. 



1.4. Non-linear effects 

To finish our summary let us briefly study possible non-linear interactions along the same line 
of reasoning. We will analyze these effects in great detail in the forthcoming sections. However, it 
is instructive to study a few simple cases which turn out to be paradigmatic examples. 



1.4--1- Shift symmetry 

Let us start by considering interactions that respect a shift symmetry for vr. Then the first 

operator we may introduce is of the form On (so that O = O). The structure of the Lagrangian 

induced by the non-linear realization of time-diffeomorphisms implies that this term comes attached 
with: 12 

- ImTTf. (45) 

For simplicity we remove the tildes from now on. Under the assumption that the linear piece 
induces dissipation (see Eq. (j20p ). it is straightforward to show the EOM becomes^^ 

VTfc + 7 - \[^^7^^i7T]k^ + a;^(fe)7rfc = -A^-^ (jOf - [d^{ddiTr)]k) , (46) 

where [ ]k stands for the convolution. We will not attempt a detailed account at this stage, but 
instead we provide some heuristic arguments to isolate the basic bits of the full computation (see 
sec. 17.11 otherwise). There are at least two effects due to the non-linearities to take into account 
(ignoring the homogenous solution, which as we discussed decays away), namely 

7r,{t) = iV-i dt'C'^it - t') I -6dl{t') - ^ dfC'^it' - t")6dUt") 

diG';{t' - i)60lii)dOlit')^ . (47) 



k 



2 



The second term in the first line comes from the quadratic term (diir)'^ in Eq. (j46p after substituting 
the forced solution for vr, i.e. vr ~ — J G60s- The piece in the second line comes instead from the 
last term on the right hand side (RHS) of Eq. (|46|) . Let us compute the contribution from the first 
non-linear term. Assuming the noise is Gaussian, e.g. 

{6df6dl5dl5di) ~ {6df6dl){6dl6di) + ..., (48) 

and using the local properties of the two-point functions together with Eq. (j37p . we obtain (for 



We will show this term arises from a —Og'^^' coupling, and we particularize to (ftvr)^ since it dominates over other 
terms, such as tt^, for fc* ~ VjH /cs 3> H. 

As we discuss momentarily, the non-linear coupling proportional to 7 can also arise from an emergent shift sym- 
metry in the non-linear response for SO. 
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fci ~ A;2 ~ A;3 ~ k)^'^ 

(7rfc7rfc7r,)(^) = J dt'dfdt'" (^^(t - t'"))' - t') (^^(t' - t"))\ .... (49) 

If we now multiply by —H'^ to transform to ^ ~ —Htt) and divide by (CC)^) using 

dtG^ ~ ~ l/(c,fe)2, (50) 
this simplified analysis indicates a value for the non-Gaussianities of order 

\M - (51) 

Unfortunately we cannot use the local approximation for the last term in Eq. (j46p since we 
assumed it applies for its time derivative, however, let us try to estimate its value by comparing 
with the one we just computed. If we take the ratio between the two at the level of the EOM we 
get (schematically) 

which supports the value of /nl in Eq. (j5ip also for this operator. 

Notice that we can also estimate the size of the non-Gaussianities by taking the ratio 



o{diTi) 



On 



clH^ clH 



C^I/nlI-^, (53) 



which is consistent with the more detailed result of Eq. (|145p . 

From here we conclude that a large value for 7 is linked to large non-Gaussianities, provided 
the operator responsible for dissipation also induces terms such as in Eq. ([15|) . or more generally 
a 7(9j7r)^ piece in the EOM as in Eq. (j46p . As we shall see throughout the paper, this is indeed 
the case for a vast class of models. 



1.4.2. Approximate shift symmetry 

Let us assume the shift symmetry vr — t- vr -|- c is softly broken by a parameter e <C 1, as it 
happens due to the slow-roll approximation. Without this invariance we can have a coupling of 
the form 

- f{t)507:, (54) 

responsible for a (local) dissipative term, plus a source noise of the form —fdOs- (As we shall see 
these terms arise from a —f{t + ■7r)0 coupling in the effective action.) Notice at linear level this 

To simplify the notation, here and elsewhere in this section we omit the momentum conserving delta functions. 
For the estimate in (I52p we used the linear part of Eq. (|47p . integrated by parts the time derivative, and used 
JdtG'; ~ 1/7. 
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operator is of the same type we studied previously, except for the overah factor of /, which we 
assume is (approximately) constant to preserve the shift symmetry. In this scenario, and assuming 
the noise is Gaussian, contributions to the three point function will be induced by 60f{t)Tr'^ at 
linear order in the response. Then the level of non-Gaussianity can be estimated to be (see sec. 
17.21 for more details) 

a^_J^Lc^;,,._m.O(.), (55) 
that is in practice very small, provided e ^ 1. 



1.4--3. Non-linear response 

Let us continue with the coupling f{t + Tr)0 but include now the response beyond linear theory, 
in which case we do not necessarily have the constraint of Eq. (|55p. Hence we have to include 
contributions to 6^"^^ Or at second order in vr which arise from the intrinsic three-point function of 
the (5Cs, i.e. {[60{z), [dO{y),60{x)]]). In general, making use of the local approximation, we have 
(schematically) 

fS^^^O§{u;)^N,goik,oj)7Tl (56) 

with goik, u}) depending on the specific dynamics of the model. From here we get non-Gaussianities 
of order 

go{k,u})TTl go{k^,uj^) 

^|^~/nlC^/nl~^^^^, (57) 

which is not suppressed by factors of f /(Hf). Unfortunately, it is not possible in general to relate 
the level of non-Gaussianities and the dissipative coefficient 7, unless the different terms in 6^"^^ Or 
are somehow related. However, there are specific situations where this happens, in which case we 
expect a connection between dissipation and non- linear interactions. 

For instance let us consider the case in which inflation is driven by a scalar field (p and the 
dynamics of the interaction is such that in the background 

O = F{$). (58) 

Intuitively Eq. (I58p follows from some basic requirement of a velocity dependence to induce dissi- 
pation. 

Now we perturb (p ^ (j) + 6(f). Given that O is a scalar, then (provided 6(p is a smooth pertur- 
bation) 



{O) ~ F 



-dcpf 



f60n ^ iVc7 ( ^ + f - li^^7r)^ + ■■■], (59) 



where vr = S(p/(j), and the factor of / appears in order to properly normalize the coupling to vr. In 
this expression we also assumed the linear piece is responsible for dissipation. (The coefficient a is 
an order one number which vanishes for the special case F{x) = |x|.)^^ One might wonder whether 



The argument applies to a generic function F{(t)), in which case 7 = 7((?1>). See sec. 17.31 for more details. 
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there is any way to get c^'s other than through y — {d(j))^ . Certainly the background breaks time 
difFeomorphisms, and therefore we have a natural timelike vector ~ d^t at our disposition. 
However, if the response of the O's is predominately determined by the field (j), we have ~ dfi4>, 
then 

dt(t> ^ n^'^^(t) = ^/-{d(t>?- (60) 

The extra terms in Eq. (j59p thus appear from the fact that the equal time surfaces set by the 
inflaton also fluctuate. This case is now similar to the one we discussed in sec. 11.4.11 and we end 
up with a 7(5j-7r)^ term in the EOM. Hence we get go — 7^^, and plugging it back into Eq. (j57p 
we obtain 

I/nlI ~ ^. (61) 

as in Eq. (jSip . Therefore, in this example large dissipation is also connected with an enhancement 
of non-linear effects. 

There is a subtle point in the above argument. As we mentioned, in general large non- 
Gaussianities do not necessarily follow from a dissipative term. We obtain large effects in cases 
where the O operators are sensitive only to fluctuations of the clock that controls the end of in- 
flation, namely the field (j) in the above example. We refer to this as having a preferred clock. We 
discuss this in more detail in sec. 17.31 

Note also that these estimations apply under the assumption of locality (in which case we have 
a well defined derivative expansion). Non-local effects can potentially increase even more the level 
of non-Gaussianities, such as it happens in the model analyzed in [27]. However in this regime the 
EFT treatment becomes more difficult. We do not explore this scenario in this paper. 



1.4--4- Non-Gaussian noise 



Going over the possible sources of non-linearties, we are finally led to consider the case in which 
correlation functions of the noise are themselves not Gaussian. If, for simplicity, we assume that 
the three-point function is local, i.e. 

{SOsii)50s{t')60s{t")) ~ iyo^6{i - t')6{t' - t"), (62) 

we get from the —f{t)50-K interaction 



(vrfcvrfc7rfe)(^) ~ j dt' {G^^{t - t'))\ (63) 



Then using Eq. ()50p we can estimate 

/N. ~ (64) 

which depends on various parameters, although clearly it can also be large. See sec. 17.41 for more 
details. 



Adding the expansion of the universe changes things a little bit, in particular we will have to 
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deal with exponential dilution. However, once we assume the dissipative effects are taking place at 
a faster pace than the Hubble expansion our results in flat space are a good guidance to understand 
the basic features of the full computation. As we shall see, most of our previous analysis remains 
essentially unchanged (provided the dissipative mechanism acts periodically over the inflationary 
epoch). 

The new ingredient is the construction of an EFT formalism from which we will obtain the type 
of terms we discussed, and more importantly the non trivial connections between the linear and 
non- linear effects. We start with the EFT setup next. 



2. EFFECTIVE FIELD THEORY SETUP 



As shown in [13!], for single clock inflation the action in the unitary gauge is given by 
S = ^ j d^x^gR+-j d'x^{p-p-{p + p)gOO) 



1 



1 



00\ 



+ i 1 d^x^gMlm + g'^'f-l j d^x./^lg^^d^^bg^'Tdp^g'''' 



- d'xV^M:{t){6Kp\l + g'^) 

where the spatially flat FRW background metric is given by 

ds^ = g^j.udx'^dx" = -dt^ + a^{t)5ijdx'dx^ , 
and the unit vector perpendicular to surfaces of constant time t. 



takes the form = —S'^{—g^^)~^^'^. The extrinsic curvature of the surfaces is 



(65) 



(66) 



(67) 



(68) 



where g^p = gj_ip + n^np is the induced spatial metric. Thus 5Ky = Ky —Hg'i is the variation of the 
extrinsic curvature of constant time surfaces with respect to the unperturbed FRW. The ellipses 
in (jGSp account for any additional term that respect (time dependent) spatial diffeomorphisms. 
Defining 



2 5S 



Einstein equations imply (a bar over any quantity denotes its unperturbed value) 



(69) 



(70a) 
(70b) 
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where H = a/a, = {8ttGn)~^- 

To introduce the vr field in the EFT we fohow Stiickelberg's trick, 

t^i = t-Tr, ^ = x\ (71) 

so that g^^ can be written as 

= + + 2d,7rg^\i){l + tt) + g^^(i)9,7r5,7r, (72) 



and 5K'j is 



+ o2fl59« - (y«s»° - H^) S„, (73) 

to Hnear order in the perturbations. Prom now on tildes will be omitted. In addition, we can 
choose coordinates so that the metric in the perturbations is given by 

ds^ = -N^dt"^ + h,j {dx' + N'dt) {dx^ + NUt) . (74) 

In this paper we will ignore tensor perturbations. 

In general the metric perturbations 5N and Ni are determined by the momentum and Hamil- 
tonian constraints. The action for tt is obtained after introducing their solution back into the 



action. For single field infiation, it was shown in 13[ that in certain regimes the metric fluctua- 



tions can be ignored, and indeed these are suppressed either by slow roll parameters, or by ratios 
of H'^/Mp. The same occurs when we include ADOF. We discuss the details of this decoupling 
limit in appendix [Pl 

The quadratic contribution to the action for it can thus be written as 



S. = ll d'xa'!^ip + p + iM^)n'-ip + p + HMl)^^-iMl + Ml^^ (75) 

(To arrive at this result we have performed integrations by parts. ^'^) The above expression can be 
re-arranged as follows (ignoring M2 and M3) 



S^= d xa —<ir - c^—^ \ , (76) 

where 



For more details on the EFT formalism see , fl^ . [2^ . As we show next, the introduction of 



The differences with are due to the fact that here we are using a term Sg'^'SKji instead of SNSEji, where 
KH = ^f^E>^. The relations between the coefficients are M? = d^M^ /2, M| = ~ S/MiHM^, M2 + M3 = 

M^{d2+d3). 
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ADOF changes the relations in (|77p . since in general they can have a non- vanishing background 
stress energy tensor. -"^^ 

3. ADDING NEW DEGREES OF FREEDOM 

In order to include dissipative effects in our system we will follow the procedure of sec. 11.21 
and introduce a set of (composite) operators that behave as an effective environment. Since we 
are dealing with gravity, we should take into account that the stress tensor corresponding to these 
new degrees of freedom, T^^ , may contribute significantly to the background. This is the case, for 



example, in the trapped inflation model of [27| where particles are created while the inflaton slow 
rolls (see sec. [8|). 

The fact that there is more than one field whose stress energy density takes an expectation value 
slightly complicates the construction of the EFT. The basic idea of the EFT of inflation is rooted 
in the necessity of having an end point for the accelerated expansion, and that there is a physical 
clock that defines a special time-slicing where the clock is taken to be uniform. This is the so- 
called unitary gauge. Time translations are spontaneously broken during inflation by the presence 
of this preferred clock, which means that there is a Goldstone boson that non-linearly realizes the 
symmetry. In the case where we add the ADOF, there is an ambiguity in the definition of the 
clock field, as the additional fields may have non-vanishing background expectation value that also 
break time-translation invariance. Nevertheless, there are two natural definitions of the field that 
interpolates for the Goldstone boson, both equally good. The first one follows the approach of \2^ . 
in which one introduces the Goldstone boson of time-translations, fr , such that the action takes 



the form ^3|, [2 



j {-Mlim'^it + + H{t + 7f)) + MlH{t + ^)5°°(vf)) + . . . (78) 

(The ellipses include other (non-Goldstone) combinations that depend on the ADOF.) In this 
approach the only tadoples (namely terms that are linear in 5g^^) are the ones associated with 
the operators \/—gg^^ and \/—g^ whose coefficients are uniquely fixed by the geometry as shown 



m 



10|, 



A second alternative, which is the one we take in this paper, is to define a unitary gauge in which 
the physical clock that controls the end of inflation is taken to be uniform. Then by performing 
a time diffeomorphism we introduce a different Stiickelberg field, that we will denote as vr. The 
main difference between the two gauges relies on the fact that now the coefficients for the tadpole 
operators, \/—g(p^ and \f--g., are not determined by the geometry and will in general depend on 
contributions from the ADOF in the background (see Eqs. (|79p ~(j82 p below). This is the case 
because we also need to include tadpole operators induced by the ADOF. The two different vr's 
are related by a mixing that involves the ADOF fiuctuations, schematically: vr ~ vr + bO. The 
field vf has a simpler Lagrangian, because the coefficients of the two tadpole terms are fixed, as 
shown in [l^]. However, it is not convenient for us because vr is not sufficient to determine the 
end of inflation. If we were working with vr, then the curvature perturbation would be related to 
the latter by a relationship of the form C, ~ -ffvr + bO. Instead, by taking the second choice, we 
have a slightly more complicated Lagrangian, yet the link between and vr is simply C, — —Htt (at 
linear order), with no dependence on the O operators. This is the case because, as we emphasized 



As we will explain later, this in not in contradiction with the results of where it was shown that the tadopole 
coefficients are uniquely fixed by H and H. This is due to a different choice for the field tt. 
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in the introduction, we work under the assumption that the ADOF do not contribute to density 
fluctuations at late times. 



This is the main difference between the cases we are studying here and the analysis of 26|, 



where additional light fields were included to the EFT of inflation of [13|]. It is easy to convince 
oneself that the two gauge choices are equivalent. We take the second. 

Let us continue with the construction of the effective Lagrangian. Following [13] our starting 
point is an effective action in a unitary gauge in which we write 



/ 



5(A(t)-c(t)gOO) + 5o, (79) 



where c(t),A(t) are certain tadpole coefficients soon to be fixed by enforcing Einstein equations. 
Note we added So = f d!^x^J —gCo to account for the dynamics of O independent of vr. (We will 
incorporate the couplings between vr and O in the next section.) Also, by construction. So is a 
scalar under diffeomorphisms. Then with 



= diag (po , a^Po > a^Po > ^^Po ) > (80) 



it is straightforward to show 



A(t) + \{Po- Po) = -M^iSH^ + H) (81) 

cit) + ^{po+Po) = -M^H. (82) 

As explained in [l^] we introduce vr following the Stiickelberg trick (see Eqs. ([7T]) and ([72]) ). 
which for the action in Eq. (j79p means that only c{t) contributes to the normalization of the 
quadratic Lagrangian in vr, so that^^ 

ciNc = 2c{t) = -2M^H - {po + Po)- (83) 

(Note this coefficient would be given by —2MpH in the absence of ADOF, in which case both 
definitions of vr would agree.) 

As an example, let us consider once again inflation described by a scalar field with the following 
action 



Stot = 1^9 (^-^idcpf - VicP) + Co 



(84) 



Since we assume the ADOF do not contribute significantly to density fluctuations at late time, our 
unitary gauge is the one where 5(p = 0. Then we obtain 

Stot = JV9 (-^'^'5°° - Vi^) + £0) , (85) 

where (^(t) is the background value. On the other hand, Friedmann equations (including the 
ADOF) tell us -$'^/2 = -Nc/2 and V{4>) = A(t), with (A,iVc) defined in Eqs. §3) and ^ (for 



Terms linear in n cancel out once the background EOM are satisfied. 
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Cg = 1). Hence the action takes the form of the expression in Eq. (|79p with the aforementioned 
coefficients. Moreover, we also get the usual Lagrangian for tt normalized by Nc, after identifying 
IT = 64>/(j). We notice in passing that assuming the stress energy tensor that follows from Cq obeys 
the null energy condition, i.e. po + Po ^ 0, then 

c^iVc < -2M^H. (86) 

Adding higher dimensional operators will shift the normalization of vr, like in Eq. ()77p . In 
particular we will generate a non-zero correction to the speed of sound, so that Cc, < 1. Therefore, 
before including interactions with the ADOF, at quadratic order the action is given by Eq. (j76p . 
with [Nc,Cs) some matching coefficients, defined as in Eq. (f77|) . 

We will not adopt any particular model for the ADOF, rather we will attempt to produce 
correlations between different observables, such as the power spectrum and non-Gaussianities, 
under some mild assumptions about the n-point functions of the type (C . . O). But first let us 
start by constraining the type of operators that we may add to the eff'ective action in the unitary 
gauge. 

4. THE INTERACTION TERMS IN UNITARY GAUGE 

We move now to the description of the type of operators that we can add to our Lagrangian 
in the unitary gauge that will induce couplings between the ADOF and the fluctuations of the 
clock. In general, the operators will have some tensorial transformation properties under space- 
time diffeomorphisms, and so they will be classified according to their rank. As it was shown in the 
analysis of [13, 26], one can write down operators containing only free upper indices. In our case, 
however, there is a subtlety we need to address since the O's are composite operators that may also 
contain the metric. Since the metric can be used to contract tensors made out of several different 
fields, we define tensor operators Oap... always with indices down, and so that 80ai3.../5g^" = 0. 

Let us give an example. Let us consider two operators, Oi = ■0^ and O2 = g'^'^d^if^dyilj, 
with '(/' a scalar field. These are both scalar operators, however, according to our prescription we 
should write: O2 = g^'^02f_iu, with 02/^i/ = d^ipd^ip. In this way the ambiguity with respect to 
metric factors is removed. Operators are then classified as a Taylor expansion in fiuctuations and 
derivatives as usual. We now proceed to illustrate the leading ones. 

4.1. Scalars 

In an expansion in metric fluctuations and derivatives, the most relevant operator is given by 

Sf = - J d^x^ hit)Ou (87) 

where Oi is a scalar under full space-time diffeomorphisms. The next type of operators can be 
organized as follows 

S? = -J d^x^g {/2(t)<5<7°°02 + hmg'^'fO^ + f4{t){6g^y04 + ...}, (88) 

where Oa, a = 1, 2 . . ., are also scalars and the ellipses include pieces involving higher power of the 
fluctuations as well as higher derivative terms such as d^Sg^^ or SK. In appendix |E] we discuss 
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briefly operators of the form 



(i^xyj—g s{t) 



(89) 



As we discussed, the operators 0x^2 may have (time dependent) background values, e.g. Oi^2 = 
Oi,2(i) + ^Oi^2-, which could also contribute to the background Einstein equations. These lead to 
corrections to T^y from the interaction between the ADOF and the one responsible for inflation. 
For example, let us take once again the example of a slowly rolling scalar inflaton and add the 
coupling 



v2 



X 



d Xy/—g 



(90) 



to a scalar field x (say we have a shift symmetry 
new contribution to T^j^ is given by 



with = p^^ = Einstein equations require 



+ c to prevent other couplings). Then the 

(91) 



-c{t) 



1 1 
+ -{po +Po) + -^{Px'P+Pxi 



(92) 



and similarly for A(t). Notice that the term in Eq. (j90p now contributes to the quadratic action 
for TT, and we get 



1 



if)-i;i.Px4>+Px<l>) 



9 



,00 



thus the canonical normalization coefficient becomes 

N, = -2MlH -{po+po)■ 



(94) 



This is a general feature: the normalization of the vr Lagrangian will be given by the difference 
between the total 'kinetic term', {p + p)tot = —2MpH and (only) the contribution from Co (the 
Lagrangian for O independent of vr). In this example Nc is not equal to 0^ but rather 



Nc = P + 2f2{t)d^, 



(95) 



where f2{t) = 0^/2 and = x^/A^. But this is what we expect upon noticing that adding the 
term in (j90p to the usual Lagrangian, —^(dcf))'^ — V{(l)), renormalizes the kinetic part of the action 
by a factor 



and therefore we obtain (using tt ~ 



1 + — 

A2 

X 



(96) 



+ 2/2(t)0,)(^ + 



TT 



ciN, ( — + 



(97) 
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as in Eq. (|96p . Note that from an effective field theory point of view, the consistency of this 
particular example with a slow rolling scalar field requires to have, even in a situation where 
^ 1 1^1, ^ s-iid therefore c^A'c remains essentially given by (p' . This is so as otherwise we 
should consider an infinite amount of terms. However this situation is not necessarily the case for 
all possible UV models, since we could instead of the term in (|9Up write a generic expansion (see 



for example 42l |) 

(xVA^) W)' ^ M4F(x2)P(X), (98) 

with X = x/^x ^^'^ ~ —{d(l))'^/M^ with M some mass scale, such that f2{t)0-^ > (j)^ . One of 
the most useful aspects of the EFT of inflation is that we do not need to worry about a specific 
realization of the background while studying its perturbations. As a result the scale Nc may be 
dominated by /2(i)0^, rather than 0^. 



Let us now return to the operator in Eq. ()87p . In principle it can also have a background value, 
i.e. ^J—gfi{t)Oi{t). Since it only contributes a piece proportional to ^/--g it can be absorbed into 
A(t) in Eq. ()79p . to ensure the background satisfies Einstein equations. Notice that the full Oi is a 
scalar so that the coupling fi{t)Oi only develops vr's from fi{t) — )• /i(t + 7r). However, if we absorb 
fi{t)Oi into A(t) then its value gets fixed as in Eq. (jSip . which we now have to expand in f + vr. 
Hence somehow the pieces from Oi{t + vr) must cancel out, and they do once we realize 50i{t) is 
not invariant under time reparameterizations and their background values are re-introduced from 

50i ^ 60i - bi{t)Ti + . . . . (99) 

In other words, the fact that the O operators have background expectation values means that if 
we split them into background plus fiuctuations, the latter shift under a time diffeomorphism.^'' 
Of course if we do not split the operator in this manner, then since O is a scalar, no tt field will 
be associated with it once we perform a time diffeomorphism. 



Let us consider for instance a coupling (jP'x? between the infiaton and a second scalar field (this 
will reappear later on), and allow for a non-zero expectation value x?{^) 7^ 0. In our unitary gauge 
we have a term in the action (f){t)'^ {x^) {t) (plus perturbations in the x's), which we can think of 
as being included in ¥{(()) (with time dependent coefficients). This corresponds to a f{t)0 type of 
coupling. As we mentioned above, we do not want to stream tt off the time dependence in x^{t) 
(because the full operator is a scalar), but this will happen once we solve for V{(t)), i.e. A(t) as in Eq. 
(jSip . However, it is easy to see these extra terms cancel out against the ones induced from Eq. (j99p . 



Let us finally briefiy comment on the slow roll approximation, since in principle the coupling 
(/i(i)Oi(t)) may break it. In general the slow roll condition can be satisfied provided 



cf [h{t)Oi{t)) ^ H 



(100) 



This requires, in addition to e = —H/H^ <^ 1 and r] = <^ 1, that any explicit function of 
time f{t) in the action, plus all background quantities associated to the ADOF, change very little 



If instead we had chosen to work with n, such that by construction aU the information about O is aheady 
incorporated in Eq. (|78p . we would still have these background values appearing in the Lagrangian from the shift 
of SO after re-inserting tt, similarly to Eq. ((99}. 
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in a Hubble time. Schematically we write: e/ = <^ 1 and = ^ 1- In practice we 

assume all the terms proportional to Oi^2it), or in general stemming from Co, are included in 
the background geometry {H,H) or into the coefficients {cs,Nc), and furthermore with their time 
dependence suppressed by slow roll parameters unless otherwise noted. 



4.2. Vectors 



Moving into vector couplings, the one with the least number of metric fluctuations has the form 

j d^x^ hit)0^5g^^, (101) 

where we have been careful in defining the vector with the index lowered as we stressed at the 
beginning of the section. As we will see when we reinsert the vr field, something unusual about this 
operator is that it only entails terms linear in vr, provided fi{t) is a constant. At higher order the 
generalization is straightforward: 

J d'x^g (/2(t)0„55™<55"° + ...)• (102) 



4.3. Tensors 



We can move on by considering generic tensors, with their indices contracted with g^^^s, as for 
instance 

I d^x^g f{t)0,...Jg'^'> . . . 6g''> = J d^x^g f{t)0'-'>. (103) 

Another type of terms, perhaps more interesting, are those coupled to the extrinsic curvature. For 
example, 

J d^x^g s{t)^d^,5K''>^, (104) 

where the factor of Mk, which we take to be much bigger than Hubble, accounts for the mass 
dimensions of Ky . We can clearly continue adding factors of 5g^^ and SKli . 



5. THE INTERACTION TERMS FOR TT 



As we already pointed out, in this paper we ignore the mixing with gravity and work in the 
decoupling limit (see appendix [D]). Therefore, following our previous sketching of the procedure, 
to construct the interacting part of the action between the ADOF and the vr's we simply replace 
(see sec. [2]) 

^00^ _l_2^_v^2^ J-(d.vr)2 (105) 
gO^^^ -5^(l + 7^) + 5flai7r. (106) 
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Also terms from the extrinsic curvature, that at hnear order induces 

5Kij a^HSij^ - didjTT + .... (107) 

As we anticipated, at quadratic order the Lagrangian for the vr field takes the form in Eq. (j76p . 

Next we include the interaction terms between the Cs and vr in the effective action. Let us 
start at quadratic level in the fluctuations. There are many operators that contribution at linear 
order in vr. However, from Eqs. <\105\ \106\ \W7\i we note that all the terms at leading order in 
derivatives can be re-grouped basically as 

AWOvr, f2it)dir, Mt)0'di7T, (108) 

where the dots include higher derivative terms. Since the non-trivial features of the non- linear 
realization of time-diffeomorphisms comes from the connection between terms with different powers 
of vr, at linear level we obtain basically all the terms allowed by rotational invariance. The first term 
appears after expanding fi(t + 7r)(D in powers of vr to first order, whereas the second term comes 
from f2{t)d5g^^. Note that there is a contribution from Eq. ([89]) to the second term, however, it 
is suppressed by a factor of H/Mk ^ 1- Terms like diTrO"^, which follow from Ofj,6g^^^, may also 
generate contributions to Cs as well as fc-dependent friction. 

For the purpose of understanding the generation of friction, we can concentrate on the linear 
order. We can therefore simply use integratation by parts and study an effective operator of the 
form 

- I <fx^ f{t)0{x)7r{x), (109) 

where f{t) provides an overall normalization scale which we assume remains constant protected by 
an approximate shift symmetry, but see sees. 11.4.21 and 17.21 In most of the expressions below we 
assume / is absorbed into O unless otherwise noted. Here O accounts for a series of contributions, 
including diO^ etc, so that we expect its Green's function to be quite generic. 

We wish to understand now under which circumstances we recover an equation equivalent to 
([2]). The main difference, as we just mentioned, lies in the expansion of the universe. A crucial 
simplification will come from our assumption of a faster than a time scale for dissipation, and 
therefore our analysis from sec. [1] remains essentially unaltered. 



5.1. Modified dynamics &i local approximations 

In sec. 11.21 we started our discussion of the effect of terms such as Ovr in the dynamics of vr 
in a flat background. The main difference now is the explicit time dependence introduced by the 
scale factor. For that reason, instead of working in frequency space, we find convenient to work in 
mixed Fourier space (t,k), keeping time as usual. Again we split the operator into pieces, 

0{t, x) = 0{t) + 50s{t, x) + 60R{t, x), (110) 

with 0{t) the background expectation value, 60s{t,x.) the stochastic fluctuations, and 60ji{t,x) 
the those induced by vr. In what follows we omit the background piece 0{t), which as explained 



This means that in practice we will work in the regime where terms like H50, that appear after we integrate by 
parts and hit the a's in the volume factors, are such HSO <^ SO. 
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in sec. 14. H we assume is absorbed in H, H, Nc, and its time-dependence is suppressed in the slow 
roll approximation {e,ri,ef,eo) <^ 1. 

Recall the first variation in Eq. (jllOp represents the noise, whereas the second one is the response 
to the perturbation induced by the vr field, and is given as the integral of a Green's function as in 
Eqs. (|14l lisp . Varying the action we obtain the EOM 

^k{t) + 3Hnk{t) + ^vTfc -^J dt'a^{t')Gg,{t,t',k)7Tk{t') = -l-60s{t,k), (111) 

with 

Gg,it,t',k)=i I -0^e-''^-y[6O{t,y),5Oit',O)]eit-t'). (112) 

(The overall normalization is given by N,, as in Eq. ()76p .) 

Our first approximation entails locality in space, and so we take the Green's function to be of 
the form 

r<0 (f f>\ 

^-(^'^''^^ = a3/2rt)l^/2(tO +Q(''^'/^^^' ^'''^ 

with Mo ^ k^,. (The factors of a^^/^ account for the fact that we work in co-moving coordinates.) 
In other words, there is a 'gap' in momentum space determined by the 'mean free path' lo ~ 
1/Mq <C We perform the same approximation for the correlation functions of the noise. For 

example for the two-point function we have 

{60s{t,k)60s{t',^)) ^ ^37^^y^^(2-)''^^'^(q + k). (114) 

To obtain a local approximation in time we assume that the characteristic time scale for the 
variation of the kernels, Tq^ <^ ^/H, is much smaller than the one of the sources, i.e. the vr field. 
(Notice that, at least in principle, is not necessarily related to Iq-) Then by changing the 
integration variable t' = t — t in the above EOM we can approximate 

1Tk{t - t) TTk{t) - nk{t)T + . . . , (115) 

The first term would introduce a mass for vr, as can be seen after using this approximation in 
Eq. (jllip . However, as we mentioned in sec. [H in this paper we concentrate in models where is 
not affected by the ADOF after horizon exit. This requires that the equations for the C modes do 
not have a mass term.^^ This imposes the condition that for a constant C the response from the 
Cs should vanish. More specifically we impose — )■ as /c/a — )• 0, such that we do not generate 
a mass term for vr. We enforce this at the level of the Green's function, imposing an emergent shift 
symmetry such that the effect of the first term vanishes, i.e. 

^ ^-Gg^{t,t-T)dT = 0. (116) 



z3/2(t) 



This is guaranteed if we have a (softly broken) shift symmetry. 
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On the other hand, the second term of Eq. pi5|) produces our desired result, where the friction 
part is given by (see Eqs. (|llip and (jll3|) ) 

^cl^-j ^-^^^^ ■ r ■ GUt, t - T)dT. (117) 

In the flat space limit this corresponds to the condition 

GZt{t, t') ^ -^N,dtd{t -t') + ..., (118) 
or Eq. (I19p in Fourier space, after re-introducing the factors of A'^c- 



The noise part will affect vr only through integrals of the Green's function whose variation time 
scale is assumed to be much longer than the characteristic scale corresponding to the noise. This 
means that vr will be sensitive only to the integral in cosmic time of I'oii, t'), and therefore we can 
approximate 

9o{t,t')c^uoS{t-t'). (119) 

If for example we would assume O is in thermal equilibrium, at high temperature T we could 
use the FD theorem which relates 

(120) 



Notice that the expansion of the universe, i.e. the factors of e~~'^, helps to improve the 
locality of the expression in Eq. (jll7p . That is to say, there is no significant influence between 
different Hubble times. In this paper we thus take 7, vq to be essentially constant up to slow roll 
effects, i.e. (^,^)~0(e). 

At the end of the day the EOM becomes 

^k{t) + {3H + 7)^fc(t) + ^TTk = -^SOsit, k), (121) 
plus the behavior of the noise dictated by Eqs. (jll4l 1119"]) . 



5.2. The homogenous solution 

Here we show that the homogenous part of Eq. (|12ip becomes negligible at horizon crossing for 
'y ^ H, which is the domain we are interested in this paper. To solve the equation we first make 
the change of variables vr = z^l'^ip^ with A = 2 + ^ 1. The equation for the perturbation reads 

(^ + l-^(l + A/2))v.. = 0, (122) 

with z = —kcsf], and r] the conformal time. Naively, as in the case with A = 2, it appears as if the 
mode freezes out when A^/z^ ~ 1, namely 2^ ~ or Cgki, ~ 7. However, as we mentioned 

already, this expectation is incorrect. In fact, we can solve Eq. (jl22p exactly and the solution for 
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vr looks like 

TTk{z) = A\y^{z) + A'^^y^iz), (123) 

where 

yi{z, u) = z'^Mz); y2{z, v) = z'^y.(z), (124) 

2^, and Jy^Yy are Bessel functions. By studying the asymptotic behavior we notice only 
2/2 tends to a finite value as z — )• 0, 

y2(^^0,zy)^-— rM~-2^/— e^P°s'^-i) for i/ > 1, (125) 
TT V vr 

using Stirling's approximation. In order to estimate 7r(z — )• 0), and consequently its contribution 
to the two-point function, we need to specify the initial conditions to extract the value of A\. The 
most conservative approach is to assume that at some given zq the mode is in the Bunch-Davies 
vacuum. (More precisely: (7r(2;o)7r(2;o)) ~ {t^{zo)'k{zq))^d.) This requires^^ 

A\ ~ for > 1. (126) 

The origin of the early time scale zq can be understood by taking our dissipative system to be 
characterized by a typically high energy scale, so that it decouples from fluctuations above this 
threshold. This implies that we can put vr in the Bunch-Davies vacuum above some this scale. 
Then as it approaches freeze out we have 

y2(z ^ 0, 1/) ~ 2^^^^ ''V2^e'^('°s^-i) ^ V2^e~\ (127) 

which is exponentially small for zq > 2i/ S> 1. Notice that even performing the matching at 
zo = 2i/ ~ A » 1, as 7/ — )• the suppression is still exponential. In other words, the mode does 
not immediately freeze out, as naively suggested by Eq. (jl22p . in fact it continues decreasing as it 
approaches z — )• y^.^^ 

Hence we conclude that the homogenous part effectively becomes unimportant for j ^ H which 
means that the two-point function can be easily dominated by the noise, as we will assume from 
now on. 



■^^ Note that the solutions in Eq. (|124l) tend to coszo and sinzo, and therefore both are required to match into the 
Bunch-Davies vaccum, i.e. e'^". This is slightly different than the analysis in sec. 11.3.11 however, notice that the 

■y\tol 

factor of Zq ^ resembles the exponential suppression e 2 in Eq. (|29[) . 
To show this more explicitly we can take the ratio ^^^y^ which, using 

'^(^'^-) ^ ^^y_^_^^ _ jrVH (f^^ ^ ^ 0, » 1), (128) 



dz 

we see goes like 



d-K jdz 



(129) 



■K I Z V 

This suggests the solution starts to deviate from the asymptotic value near z ~ y^, or Csfc* ~ \/lH , as we 
anticipated in sec. 11.3.11 
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6. THE POWER SPECTRUM 

The computation of the power spectrum fohows the same step as in sec. 11.3.21 except that we 
have to deal with a somewhat more elaborate Green's function. The particular solution of Eq. 
(jl2ip is given by 

kcs n 60s 
'^kin) = irrjrl I drj' G.i{kcs\r]\,kcs\ri'\)— (130) 

where 



G^{z,z') = ' [Y,{z)Mz') - Mz)Y,{z')] , (131) 



2 \z' 

with z = —kcgT] and z' = —kcgr]'. Then with the use of Eq. (I114| [TT9]l we obtain 

P^{k) = {TTkT^k)o = Y dz'{G^{z, z)f. 

For kcsT] —7- and kc^riQ — ^ — oo we find 

4 



i6*(^ + i)3r(^) ^ ^i^hVHJ^- 

" N^ikcsf 
or (P^ = H^P^) 



^^^^^^ W^^F(fTif^ " ^'W^ ^'''^ 



Hi 



A, = k^P,{k) ^ -hV^^^^-^^- (133) 



(Recall the * means that the quantity is evaluated at freeze out Csk/a{t^) ~ VT*^*-) 

In Fig. [1] the power spectrum is shown as a function of 7 and also |A;Csr/o|- Notice that the 

dependence on rjQ drops out once we take \kcs'qo\ ^ 1-^^ 

If we were to assume the ADOF are in thermal equilibrium at a (high) temperature T, using 

the FD theorem we would obtain 

rp tt2 

k^CkCkh ^ xA^ ^ , , J.. . ■ (134) 
2c* {c*/Nc) 

These are exactly (up to numerical factors) the results in Eq. (j^ [33]). We note that an 



expression similar to Eq. (jl34p was first introduced in 3J, |35|], and plays a key role in warm 
inflation models [36]. 

A few comments are in order. First of all, since we lump a series of operators into a single one, 
i.e. On, the local approximations look significantly different depending on the type of terms we 
are dealing with. However, once our assumptions are enforced the analysis is quite robust, and 
does not depend on the very details of the models. The only constraint we need to ensure is the 
scale invariance of the two-point function, which is guaranteed as long as 7, i^o remain relatively 
constant during inflation (more below). 



Note that even in the case where 7 <C this contribution is still larger than the result for the Bunch-Davies 
vacuum, provided i/o > ^NcH^. Hence, as long as A'^c < eH^AIp, it is enough to have vo eH'^Mp. 
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1.0 h 




5 10 15 20 

FIG. 1: Top: Dependence on 7 of power spectrum given in Eq. (|132[) normalized as p-ni'y) = 
QN^ {kcs)^ PTr{k) / {tti/o) , for zq — > +00; Bottom: The power spectrum Pnik) for \kcs'q\ — as a function of 
zo, for = (solid), = 3 (Dashed), and -f/H = 9 (Dotted). 



The distinction between models starts to play a role when we move to the non-linear level. For 
example, the friction term could derive exclusively from the operator O^J^^*^, which does not lead 
to non-linear couplings between O and vr; or it may be produced by the scalar couplings f{t)0 or 
OSg^^, in which case we do generate quite distinct non-linear interactions. As we discussed before, 
OSg^^ appears to be the cleanest, for which the value of /nl is tied up with the friction coefficient. 
This is not always the case for f{t)0, where the level of non-Gaussianities depends on the details 
of the model. However, as we argued before, in cases where there is only one preferred clock and 
the (non-linear) response 60r has an emergent shift symmetry and /nl is equally enhanced. We 
make this analysis more precise in what follows. 

7. NON-GAUSSIANITIES 

The non-Gaussian features of the models we discussed in this paper are perhaps the most 
interesting results from the point of view of observational signatures. Indeed, while the two point 
function is essentially fixed by the symmetries of the quasi de Sitter space, and characterized just 
by two parameters, namely the overall amplitude and the tilt, higher order correlation functions 
enjoy a functional freedom that, if observed, it would allow us to decode a much larger amount of 
information about the physical mechanisms that produces them. 

Studying the full spectrum of non-Gaussian signatures for the vast class of possible operators 
we introduced is difficult, mainly because of its generality, and thus lies beyond the scope of the 
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present paper. Therefore we will concentrate on the leading (scalar) operators 06g^^ and f{t)0. 
Very interestingly, we will show that enforcement of the symmetries allows us to relate the linear 
and (some of) the non-linear contributions in the EOM, such that large dissipation will be linked 
to large non-Gaussianities. In the next sections we assume dissipation takes on the form 771". We 
will study the case of higher derivative dissipation in sec. 17. 5i 



7.1. OSg°° 

Here we analyze a situation in which the dominant interaction term has one time derivative of 
TT, and it is produced by scalar operators of the form —06g^^. This term has a piece which is linear 
in TT, more precisely 2071". (In what follows we absorb the factor of 2 induced from Eq. (jlOSp into 
the operator O, i.e. O 0/2.) 

Notice that now we have a choice about the local approximation described in sec. 15. H namely 
it may apply to (SOsSOs) and G^^, or {SOs^Os) and G^^- I^'o'^ ^I'st case we will not obtain 
the usual dissipative term 771", but higher derivative terms (as in the ADL force). Here we analyze 
the second possibility and discuss the former in sec. 17.51 

Since we take the scale of time variation for 50$ to be much shorter than the main 

contribution to the noise comes from 60s ^ H50s- We can then basically follow the same steps 
as in sec. l5.1l to show that to linear order the equation for tt reduces to Eq. (jl2ip with 60s replaced 
by dOs- Then the power spectrum is given in Eq. ()132|) . with the replacement i^o ~^ ^q- More 
explicitly, to second order the equation for tt is given by 

^ + 3Hir - -^TT = -N^^ i 60{l + n) + 60n - -^di{80di^) \ , (135) 

where 60 = 60s + 60r. From Eq. (I105P we see that the force disturbing 60 is given hy F = 
TT + 7r2/2 — (5i7r)2/2. Using however the local approximation for the time derivative of the response 
part, 60r as in sec. 15. 1^ we obtain 

<^C)..Ar.7(^ + y-^). (136) 

As we already noted in sec. 11.41 working within this local approximation (and without making 
additional assumptions) we are not able to compute all the contributions to non-Gaussianities, 
since in Eq. p35p there appear not only 60 but also 60. In general, in the perspective of EFT we 
do not foresee any fine cancellations and, as we estimated in Eq. (j52p . we expect the contributions 
to the level of non-Gaussianities of the non-local terms to be about the same as the local ones. 
Therefore, in what follows we concentrate on the terms involving only O and reduce Eq. ()135p to 



c^V^ 7 diTidiTT 37 



7f + (3i/ + 7)^ - ^^7r = -ir^^ - -^ir' - N^'60s{l + ii). (137) 



Let us start by focusing on the contribution of the first term on the RHS: '^{di'ir)^. To analyze the 
non-Gaussianities we decompose 7r = 7ri +7r2, where the subscripts represent as usual the order of 
the solution for the fluctuations. Then 



kc n 

TTi{k,rj) = jy ^2 I dr]'gj{kcs\'n\,kcs\ri'\)60s, (138) 
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where gy{kcs\r]\, kcs\r]'\) = Gy{kcs\r]\, kcs\ri'\) / {kcgrj')'^ with defined in Eq. p3ip . and (using 
Eq. (dMl)) 

7r2(fc3,0) = 2^2^^ J dr]'r]'^g^{0,k3CsW\) J ^^;^|q||k3 - q|(q • (kg - q)) 

rv' rv' 
X / dr]" g^{qcs\r]'\,qcs\r]"\) j dr]''' g^{\\^'i - q\cs\r]'\,\\<iz - q\cs\r]'''\) 

X 6ds{ci,7]")6ds{'ks-q,7]"'). (139) 

(Notice 7r2 is sourced by terms quadratic in the perturbations: 50s, t^i-) 

We want to compute the 3-point function for C ~ —Htt in the hmit — )■ 0: 

(C(ki,0)C(k2,0)C(k3,0)) = (Ci(ki,0)Ci(k2,0)C2(k3,0)) +cychc sum in h's 

= {2Trf6^^\ki + k2 + 'k3)F{kuk2,k3). (140) 

Then, after Eqs. ()114p and (jllOp plus the condition that the noise is Gaussian (we will consider 
non-Gaussian noise later), i.e. 



{60s{k2,v')SOs{k,,fj)60s{ci, Ti")60s{k3 - q, r?"')) = C^^f^o^^^jFJ^^^ (141) 
X \^6{v' - r]")5{fi - ?/")5(^Hk2 + q) + - r]")5{i - T]"')5^^\ki + q)} (for ^3 / 0), 
defining Xi = ki/k, with k an arbitrary scale with units of momentum, we obtain (for rjQ — t- —00): 

F{xi,X2,X3) = k^F{ki,k2,k3) = -——^xlxlx3{xl- xl- xl) / dyy'^g^{0,X3y) (142) 

^^^c Cfi Jo 

+00 r+00 



X / dz z gj{x2y,X2z)g^{0,X2z) / dww g^(xiy, xii(;)5-y(0, xiif) + cyclic sum in Xj's 

(we performed a change of variables: y' = —kcsij', z = —kcsf]", w = —kcgr]"). 

The associated parameter (defined for equilateral configurations) is given by 

F(1,1,1) = ^/^1A^, (143) 
(with = k^P(^), and in terms of the Green's functions: 



•'NL ~ ^2 



(^+1)V(^ 



+00 f+00 

dyy'^g-yiO,y) / dzz^g-y{y,z)g-y{0,z) 
Jy 



X / dww^g^{y,w)g-y{0,w), (144) 



y 



where we have used the power spectrum given in Eq. (|132p (with uq — t- f »=,). Performing a numerical 
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integration we extract 



in the strong dissipative regime 7 S> i?. Notice this can be large even when Cg ~ 1. Using the 



analysis from Ref. 15| we can also quote a bound at order of magnitude level: 



^ < 500. (146) 



In Fig. [2] we plot the shape for two values of 7 (7 = AH and AQH) . Notice there is a peak on 
equilateral configurations, and another peak around X2 — 2:3 ~ 1/2. The latter is one of the main 
features in the bispectrum, and it is due to the fact that the fluctuations are dominated by the 
noise. 



For moderate values of j ^ H, the shape resembles the orthogonal one described in 15|, which 
is currently at 2a level in the WMAP 7-yr data (43I ]. 

Similarly, for the second source term on the RHS of Eq. p37p we obtain 



5 7 
2H 



/2-87//^vr2r(^+4)'\ .+00 .+00 
; -§- / dyg^{0,y) / dzz^y—^{y, z)g^{0, z) 



neq 

- ' '•■—{y,w)gy 



X / dww%^{y,w)g^{0,w). (147) 



y 



Despite appearances, the contribution from this term to does not increase with 7 and becomes 
an order one sub-dominant effect. This follows from the properties of the Green's functions, and 
in particular because time derivatives scale as powers of H (rather than yJ^H), since this is the 
time dependence of the Green function at freeze out. This is indeed what we expected judging 
from our results in flat space of sec. 11.3.21 ("Tcq" ~ ^/H). See appendix [Fl for some collective 
details on these Green's functions. 

For the last term in Eq. ()137p we find 



req 



^^^j 1 /■+°°.,.,.4.„ ,.^^2 f^^ ....A. <h 

dy 



dyy\g,{0,y)f / dz/y-p(y, z)g^(0, z), (148) 



which is also a contribution of order one that does not increase with 7. The shape corresponding 
to this term is shown in Fig. [3] for 7 = lOff. 

As expected (because the contributions involve time derivatives of tt, and also because 50 is 
local) all of these shapes are suppressed as we approach the squeezed limit. That is when one 
of the momenta, say ki = ki (long mode), is much smaller than the other two (short modes), 
i.e. k^ <^ ks (with ks = |k2 — k3|/2). However an interesting question remains, and that is 
whether the suppression entails one (or more) power(s) of k^ compared with the local shape, 
since this particular scaling could potentially distinguish signatures from single field models in 
measurements of the scale-dependent bias We return to this point in the following section. 
(Let us point out that from Fig. [3] the shape itself approaches zero in the squeezed limit ki/ks — )• 0, 
thus F(xi,X2,X3) cannot scale like 1/x^ = {ks/k^)'^, but rather as a softer power.) 
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X2 




FIG. 2: The shape F{x2,X3) = xjxj^^^^^ given by Eq. for 7 = 4i/ (top) and 7 = 40i? (bottom). 

To avoid showing equivalent configurations twice, the function is set to zero outside the region 1 — X2 < 

X3 < X2- 

7.2. f{t)0 I: linear response 

Let us now consider the contribution from an interaction of the form f{t)0 as in Eq. (j87p . We 
divide the possibiUties in two: hnear and non-Hnear response. We treat the latter in the next 
section. 

To compute non-linear effects in linear response theory we use the force affecting the ADOF to 
second order in it. This is given by F = /vr + /'7r^/2 + . . ., where the dots stand for terms with 
higher powers of vr proportional to higher temporal derivatives of f{t). As we mentioned in sec 
11.4.21 ill general we neglect these terms in the slow roll approximation. Therefore in what follows 
we will treat / as essentially constant. 
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1.0 




0.0 



FIG. 3: The shape F{x2,X3) = xjxl^^^^^ for 7 = lOH corresponding to the last term in Eq. (fW)) 

Working within the local approximation for G^^, the response part SOji can be written as 

60b^ ~ ^ ( -k + ^irn ) , (149) 




up to terms suppressed by higher powers of the slow roll parameter (s). (Recall we assume the 
presence of an emergent shift symmetry to neglect terms which do not involve derivatives.) Thus 
the equation for n to second order reads 

n+{3H + 7)7T + + 2jl Mfc = -N-^ (f60s + f[7rSOs]k) , (150) 

where the noise part 50s satisfies Eqs. ()1141 1119p . By comparing this equation with Eq. ()12ip . 
it follows that the power spectrum is the same as the one given in Eq. ()132p . up to an overall 
factor of namely vo ~^ ^fO = f^^o- Therefore, to guarantee the scale invariance of the power 
spectrum we require i'fc>/{vfc>H) ~ 0(e) (more on this below). 

As we stressed in sec. 11.4.21 the contribution to /nl from the non-linear terms above will be 
suppressed by f/{fH) ~ 0{e). Notice that at this order we can no longer ignore other effects, 
such as the mixing with gravity. Moreover, we show in the next section that non-Gaussianities will 
be dominated by the non-linear response. Nevertheless, we include here the computation for the 
terms in Eq. ()150p for two reasons. First of all we explicitly show that even though the non- linear 
terms are proportional to 7, it factors out in the final answer; but more importantly because it 
appears to give us a contribution in the squeezed limit which could be non-negligible in cases where 
/ is not a constant. 

Let us start analyzing the last source term in Eq. ()150p : fndOs- Since this term does not involve 
enough derivatives of tt, it is easy to see that the shape, given by F{x2-,x^) = a^i^i'^F^frff^' ^^^^ 
have its maximum contribution in the squeezed limit. In such limit the parameter /nl is given by 



,3^67x2^1 6 {P(^{i)Pi^{x2) + Pc(i)Pc(x3) + Pc(a;2)Pc(^3)) ■ 
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Then following a procedure similar to that of the previous sections we find 

whatever the value for 7. This agrees with our previous estimation in Eq. (j55p . Similarly, one 
can compute the contribution of the last term on the LHS of Eq. (jl50p . Even though this term is 
proportional to 7, the parameter /nl does not result in a significant increase compared to the one 
above. This follows from our heuristic arguments in sec. 11.4.21 but can also be seen directly from 
the properties of the Green's function, see appendix [Fl Since the approximate shift symmetry 
requires f/{Hf) ~ 0(e), this leads to very small /^^ '^^ (H52p . as we mentioned before. 

One may nonetheless worry about the result in Eq. (I152p had we assumed f,iyo were not 
approximately constant, while keeping constant the product i^fo = f'^vo-, which is the combination 
that appears in the power spectrum. However, notice that the above computation is incomplete 
since at non-linear level there are contributions that arise from the fact that tt also affects the 
probability density functional for the noise. (These are suppressed in the limit Oq <^ Hvq.) If all 
the contributions are taken into account, at the end of the day we should have /^^ proportional 
to {us — 1) to be consistent with the fact that we did not include any type of perturbation outside 
of the horizon other than the clock [l^ . |3] • 

To better understand the properties in the squeezed limit, let us compare the contributions 
from (JOtt and 50tt type of non-linearities, which appear commonly in our study. It is not difficult 
to see that if we fix the O operators to be the same in both cases, the relevant difference between 
the two is that while the shape for the latter contains g^{xLy, z), the former involves ydyg-y{xLy, z) 
(where z and y are integration variables). Using 

ydyg^{xLy,z) xly'^ 
g^{xLy,z) ~(l + 7/^)' ^ ' 

for yxL ^ 1, and taking into account that the integrals are dominated by values of y and z near 
y^j/H (see appendix |F]) , we conclude that the squeezed limit for SOn is suppressed by a factor of 
x'j^ = (ki/ks)'^ with resp ect to the contribution from SOtt, in the region x| <C H/'y. This agrees 
with the findings in [45(]. We will discuss the squeezed limit and consistency conditions in more 
detail elsewhere. 

We study next the non-linear response for a f{t)0 coupling, which turn out to provide the 
largest contributions to non-Gaussianities. 



7.3. f{t)0 II: non-linear response 

So far we considered the response 60r to linear order in vr. This entails the knowledge of only 
the two point function G^^- (see Eq. ()112p ). However, we can also induce non-linearities by going 
to higher n-point functions of the 60's, for instance including their three point function 

CZ{x, y, z) = {[60{z), [60{y),60{x)]])6{ty - U)e{t. - ty), (154) 
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while keeping the leading order force -F('7r) = fir (in what follows we treat / as constant). Then 
we have 

f60R{x) = -p J Gg,7riy)dy + fj dydz Cg,{x, y, z)n{y)7,{z) + ... (155) 

where the ellipses represent higher order corrections. 

A priori, even assuming both are local in space and time, we would not expect any relationship 
between and C^^- However, as we argue below, we will have a connection in cases where the 
response of the ADOF is governed by a single (preferred) clock. 

As we emphasized throughout the paper to obtain a dissipative term from an operator f{t)0 
we require (assuming a local Green's function) 

/(5«0^j = 7Ve7^. (156) 

Notice that this expression can be suggestively re-written as 

= iV,7 ^^'^ {n^d^^it + ^)} , (157) 

where (note n° = 1 up to 0(7r^)) 

ra^ = t = t + Tr, (158) 

yJ-g'^PdJdpi 

such that Eq. ()157p makes the symmetries manifest: O is a scalar operator. This allows us now to 
extrapolate the response to all orders in vr, that is^^ 

= A^,7 '^^"^ {n^d^it + tt)} = 5^") (^^J -g-PdJdpi^ . (159) 

It is then straightforward to show the above equation implies a relationship between linear and 
non-linear terms in the EOM. For example in the unitary gauge, i.e. t = t, 



fS^^^Oji = N,j5^''^^/^, (160) 

which allows us to read off the tt interactions using g"^ (vr) = -1 - 271- - 71-2 (5.7r)2. Hence 

= -iV,7^, (161) 
-M^^Or = lN,j{dinf. (162) 

We thus get a non-linear term that resembles the one we obtained for the case of an exact shift 
symmetry. (See the first term on the RHS of Eq. (jl37p . compare also with Eq. (jl36p .2'') Then the 
analysis of the non-Gaussianities follows as in sec. 17.11 (see also sec. I1.4.ip . so that 

hLC^-^, (163) 



^® In principle this can be generalized to Fo (^\/ —g"''dvtdpt^ . For simplicity we restrict to Fo{x) = x. 

'^'^ Notice that for the type of non-linear response in Eq. (|160p we did not obtain a tt'^ term. However, these can be 



generated if we allow for a generic function Fo{\/—g^)- 
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as anticipated in sec. 11.4.31 The shape will also resemble Fig. O 

The previous argument is perhaps better illustrated if we particularize to the traditional model of 
scalar field inflation. The crucial point is that the response to the perturbation must be compatible 
with diffeomorphism invariance. Then in order to induce dissipation we need factors of (p which in 
turn lead to large non-Gaussianities (for 7 ^> H) as sketched in 11.4.31 

The reasoning goes as follows. Let us consider an interaction Hamiltonian that takes the form 
(with O a scalar operator) 

i^int = 00. (164) 

The EOM for the inflaton thus becomes 

D^(t> = O, (165) 

where -D^<^ = represents the EOM in the absence of O. There is certainly a dynamical system 
behind O, with its own Co- Nevertheless, in the spirit of EFT we do not make any assumption 
other than the coupling in Eq. (|164p . We start now by solving for the background, i.e. <^(t). 
Assuming a local response we have 

{m^) = Fo{4>,h, (166) 

where the brackets emphasize we are computing the response in the background given by (p. Plug- 
ging into Eq. (I165P we get 

D^^ = Fo{^,$). (167) 
In order to have the type of (velocity dependent) dissipation we study in this paper we will impose 

Fo{^J)c^\$\. (168) 

We discuss a specific example of this behavior in sec. [8] (albeit with a more elaborate response 
function). The next step is to perturb (p ^ (j) + 6(f). The crucial piece of the argument is what 
replaces Eq. (|166p . There are several ways to attack this question. Perhaps the simplest is the 
following. Since O is a scalar operator (by construction) the expression in Eq. ()166p must be 
covariantized. (Notice Lorentz invariance is only broken by the presence of a preferred frame, i.e. 
in Eq. (llSSp with t replaced by (j).) Therefore, as we argued, we should rewrite Eq. ()166p as 

{(P\0\cP) = Fo (0, n^^d^cp = -g^^'^d.cpd^cP/^-g^'^d.cPd^cp) , (169) 

where we assumed the function Fq is not essentially modified, which is a valid (local) approximation 
provided 54> is a smooth long wavelength perturbation.^^ 

Once again we go to the unitary gauge, and choose coordinates such that <j) = (p, i.e. 5(j) = 0. 



In other words, we need the scale of variation of the extrinsic curvature of constant time surfaces (recall K^-' ~ 
V(3)n^) to be much larger than the typical length scale M^^ , which controls the local approximation for the 
interaction between (f) and O, such that we have a well defined derivative expansion. For instance if we take the 
perturbations at horizon crossing, namely 5(fn:^ , and say Mq ~ g\(j>\ with g some coupling constant, then we require 
fcj <C Mq, or 7 < g\<f)\/H. (In all generality, even if suppressed, we should in principle include also corrections 
induced by SK^^^.) See sec. IB.ll for more details. 
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This means 



,00 



(170) 



We then introduce the perturbation 6(j) (our vr) using Eq. (|167p and the Stiickelberg trick. On the 
LHS of we get the usual term, whereas on the RHS we use the requirement of Eq. (|168|) . then 



similarly to what we obtained from Eq. (jl60p . with vr = 6(f)/ (j). Then, as advertised, we get a 
dissipative term 771" and also 7(5j7r)^ with its subsequent enhancement of non-Gaussianities as 
shown in sec. I7.1[ 

There are some subtleties in the previous discussion that must be stressed out. Strictly speak- 
ing, the presence of dissipation requires only a linear coupling to vr. In our discussion the non- linear 
term obtained from n^d^ir in Eq. (jl59p was intrinsically related to n'^ being the time-like vector 
orthogonal to the uniform slices of the physical clock that controls the end of inflation. However, 
in general the presence of ADOF may also include time-like vector operators, let us collectively 
denote them as n^, that might have non-zero expectation values in the background. Hence, 
rather than n^S^vr we may have u^Sj^vr in Eq. ([EH). 29 In this case there will not be a guaranteed 
relationship between the linear term 771" (now derived from {uq) = (1,0,0,0)) and the non-linear 
counterparts. In this situation large dissipation may not necessarily lead to strong non-Gaussianity. 

It is worth pointing out that had we chosen to work directly with the Goldstone boson tt that 
we introduced in sec. [3l we would have had no vector-like operator taking a non-zero expectation 
value. This means that (up to normalization) the only vector taking a background expectation 
value is n'^ ~ g^"du{t + vf), which appears to induce the type of term in Eq. ()162p . However the 
two descriptions are equivalent, and in fact, upon noticing (schematically) tt ~ vr + iJO, even in this 
case we cannot in general prevent an interplay between vr and fluctuations of vector operators Sv!^ 
such that the term in Eq. (jl62p does not get generated at the non-linear level. Therefore, in order 
to ensure a connection between dissipation and non-Gaussianities, we will demand that there are 
no such cancellations. 

There is no simple way to implement this condition, although in practice this amounts to 
requiring that the 50^s are sensitive mostly to the field vr whose fluctuations control the end of 
inflation. In this paper we referred to this condition as having a preferred clock. 

Hence we conclude that the existence of such clock inevitably entails large non-Gaussianities in 
the strong dissipative regime, either from Og^^ , as explicitly shown in sec. 17. H or via non-linear 
response for f{t)0 as we just described. We will study relaxing this assumption in future work. 



As a final source of non-Gaussianities, we have to consider the effect of the intrinsic non- 
Gaussianity of the noise fluctuations. Let us consider, as in Eq. ([62]) . a non-Gaussian noise 50$ 



This may happen if the response of O is controlled by Uq. 




(171) 



7.4. 



Non-Gaussian noise 
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that is local, that is 



Then using 



fkcs n 

7r(t,k) = ^ I dr]' g^{kcs\v\, kcs\r]' \)60s{v' ,k) 



Vo 



we find (for r/ — )■ and r/o — )• — oo) 



F{xi,X2,X3) 



+00 



jy3^6 ^1^23:3 J dzz g^{0,xiz)g^f{0,X2z)g^{0,X3z). 



(172) 



(173) 



(174) 




FIG. 4: The shape F{x2,X3) = x^xl^^^jfl^ obtained from Eq. for 7 ~ lOH. To avoid showing 

equivalent configurations twice, the function is set to zero outside the region I — X2 < < X2. 

The shape is plotted in Fig. [H As expected there is a peak on equilateral configurations, and 



JNL 



(175) 



o 



as we anticipated in Eq. ([& 



7.5. Higher-derivative dissipation 

Let us consider now the same situation as in sec. 17. H where the dominant interaction has one 
derivative of vr: 5g^^O. However, unlike before, we will now assume a local approximation for the 
two-point functions of O (or ImG^^ ~ w), such that we get 



50 



R 



o 



(176) 
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where Tq (with units of mass) depends on the specific details of the model. Given that the 
interaction term goes like tt6(D, the EOM becomes^" 

vr + ^vr + 3/7vr - (a vr) = ^3^, (177) 

Note this equation is now reminiscent of the ADL force. The presence of higher derivative terms 
forces us to restrict w ~ <C To, and treat the higher derivate term as a small perturbation. In 
this case freeze out occurs at u) ^ H and the leading effect is not on the two-point function, that 
here is dominated by the vacuum solution, but on the three-point function. This can be estimated 
by comparing the non-linear term in the EOM to the linear term. The second order equation can 
be written as 

vf2 + ^vra + i/ f 3 + ^ j ^2 = Jfc, (178) 
where the term contains the quadratic sources. From the coupling between vr and bO we have 

= "iVe 1 a3 + /' (^^^^ 

where bO = 60s + SOr, and the subscripts represent to what order in vr the contributions are to 
be computed. For estimating purposes we consider the last one, which dominates for small speed 
of sound models. After using (jl76p the resulting level of non-Gaussianity is given by 



JnlC ~ 



vr 



H^tt'^ 1 H 1 1 , , 

ciTo To cj cj 



which tells us that this can be at most a subleading signal. 

8. MATCHING DISSIPATIVE EFFECTS DURING INFLATION 

Even though from the EFT standpoint the conditions on and noise kernels may be taken 
as a given, it is still desirable to be able to identify those situations where we can be assured these 
conditions hold within certain degree of approximation. For instance we expect such description to 
be valid in cases where the memory effects decays sufficiently fast, as in Drude's model of Eq. ([7]). 
This would be the case provided the dynamics of the Cs is such that the back reaction on the 
vr field becomes negligible after a (short) interaction time. In this section we study a specific 
realization in the spirit of the trapped inflation model of [27I], and also briefly discuss the warm 



inflation paradigm ||34l436l| . See appendix [H] for a class of models with 771" dissipation with a On 
coupling. 

8.1. Local trapped inflation: f{t)0 J^ti^l^ ~ <Pt)'^xf + ^Xi) 

We are after an example in which we can apply the local approximation in time. In trapped 
inflation the energy of the inflaton is being transferred into the ADOF both because the particles are 
created and also due to the increase of their mass with time. In the original model this is important 



We also generate corrections to the speed of sound of order H/To that we neglect here. 
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until the particles have diluted enough due to the expansion of the universe. The characteristic 
time scale for this to happen is given by H~^, which does not allow us to use the local response 
functions that we used so far. However, this characteristic time scale can be significantly reduced 
if the produced particles decay into yet other degrees of freedom, with decay rate much faster than 
Hubble. With this in mind let us consider the trapped inflation action 



'trap 



-d^^df'<p-v {<!))}, (181) 



where (p is the inflaton and V{4') is its potential, and we add an interaction term, 5int(Xj; V')) 
that characterizes the coupling between Xj's and some additional degrees of freedom (p that leads 
to a decay rate for the Xj's satisfying F^- ^> H. The only condition we need to impose is the 
requirement that S'int does not modify the leading order picture from S'trap- In particular we need 
to ensure it does not generate a large mass for Xi when (j) (pi. For example we can have 

I ^ 

with 93 some scalar field which we assume has some small mass, ni^p ^ \J g\(l>\, and does not couple 
directly to (j). In what follows we show that this model admits a local approximation. (To simplify 
notation in this section we drop the bar for the unperturbed values of (p.) 

As in the original trapped inflation scenario, particles associated to the field Xi ^^re created 
when ap proaches (pi. The computation of the production of particles follows the same line as 
in 



271 . |28| |. The time interval 5tc = t — ti during which modes of a given Xi held do not behave 



adiabatically (so that the corresponding particles are produced) can be estimated as 



dtc^l/Jg\<P\, (183) 



which we want to be shorter than a Hubble time, 6tc H ^, hence 27 1 

H^<-9\^\. (184) 
Note this condition implies that after particles are produced the mass satisfies 

M^AU + ^tc) = g\(l) - cf>i\ g\<p\6t, > H. (185) 
Moving to the decay rate, using Eq. ()182p we can estimate 

jc 
\? 



rx. ^ (186) 



To ensure the decay time scale is longer than the time it takes to create the particles we need 

1 ..imf" 



- > F^, ^ ^ k^>^g\<p\, (187) 

where we have approximated M^- ~ ^|0|r^. = M^.(r^.), with r^. the hfetime of the X2-P^i'ticles, 



44 



while at the same time F^^ ^ H, which is necessary for our local approximation.^^ After the 
particles are produced, the occupation number for each species is given by 



= e -'(h)-'ih)^ (188) 



where K{ti) = y g\4){ti)\, and the number density 

n^^iU + = M^al ^ (189) 



Then, taking into account the decay rate of the particles for later times, we have 

nx.it,U) - ^ ( 4^Ve-rx(*-*.) Q^t-U), (190) 



(2^)3 \a{t) 

where in order to obtain an estimation we have approximated F^. by a constant T^. The EOM for 
the unperturbed inflaton becomes 

4> + 3H^ + V'{<P)+gY,nx.it,t^) = ^■ (191) 

i 

As the sum over the particle production events is difficult to deal with, we will replace it by an 
integral. Defining A = (pi+i — (pi, the last term on the LHS of Eq. (jl9ip can be thus approximated 
by 

E-x.(M.)- (192) 

This is a good approximation only if the variation of the integrand is small between production 
events. This is quantified by the conditions 



i3H + T^)\A\ ^ |0A| ^ 



Hence, for T^.^ H and \(j)\ <^T^\(j)\ we can approximate the integral as 



(193) 



oTn (194) 
9 2_,n^AW- (27r)3 - r^|A|(2^)3- 

We discuss the consistency of our approximations in more detail in appendix HI 
We now define the (composite) operator f{t)0 = J2i fi(t)Oi, where^^ 

/,(t) = ^(</)(t)-0,)2, 0,=xl (195) 



Notice that r^^ > 5tc requires V^i < Mxi{Txi), which also consistent with the particle interpretation. 

We chose to work with f{t)0, rather than each individual fiOi, in order make direct contact with our analysis of 

the EFT in previous sections. 
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whose background expectation values are given by (see Eq. p89|) ) 

him = y Y.(^{t) - c^^?xl = \Y. (196) 

i 



Then, using Eq. (1194^ . the expression in Eq. (I19ip becomes (where we add the gradient piece, 
that vanishes in the background, for later convenience) 



^ + 3/f0-|a?^ + n« + f^f5;^ = O, (197) 

which is local in time as we advertised. (The same applies for the EOM of the perturbations we 
study in the next section.) 

To wrap up this section let us comment on the emergence of the shift symmetry, which plays 
an essential role in our analysis. Even though the EOM has a term ^^^^^(^ — (t>i)\4 which is 
apparently not invariant under <j) ^ 4> + c, secretly it is, since we end up with Eq. (|197p which is 
manifestly invariant (ignoring the shift in V{4>) which we assume is small). 

Notice, first of all, the non-perturbative result 

Xl - (198) 
9\(P - 9i\ 

which cancels the explicit factors oi (p — cpi. However this is not enough, since n^. depends on ti (the 
time defined as 4>{ti) = (pi), which in turn depends explicitly on (p. Nevertheless, this dependence 
is ultimately removed by the presence of the sum, which is the key feature. At this point the EOM 
becomes invariant under (p ^ (p + c. (This is the case because we can absorb the shift in (p into 
a redefinition of (pi, which is summed over a large number of periods.) We call this an emergent 
shift symmetry. 



,2 

s a2 



8.1.1. Perturbations 

In order to obtain the equation for the perturbations at linear order we expand (p ^ (p + 5(p \n 
Eq. ()197p . Taking into account the contribution of the noise we obtain 

6(P + ^5(P + 'iH8(P + V"{(P)5cp + -^^-^^L_50 = _^An^, (199) 

with 

An^ = 5^M^,Ax,^ (200) 
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the variance of the number of particles produced. Hence we arrive at an equation as in with 
an effective damping rate 7 given by 



5 g5/2|^|3/2 

2r^|A|(2^)= 



7= o T^ IAI.n-N3 - (201) 



In this expression we ignored the spatial variation of the field in An^^ at linear order, provided 
the physical wave vector /cph = k/a satisfies k^, <^ k, and also fcphA <^ in order to be able to 
replace sum into integrals. 

Notice that while each individual /j is not a constant the power spectrum still obeys scale 
invariance. This can be seen by redefining Oi — )• n^^,., in which case the new fi{t) ~ gMy..{(p) does 

obey fi/ifiH) <C 1 in the slow roll approximation. Indeed, the response at linear order in 5(f) 
becomes 

e. 5/2|i|3/2 

where we used 5(f)/ ^/Wl — and defined / = g(f>, with the normalization Nc ~ as explained in 
sec. [3l Likewise, for the noise (see Eq. (|200p ) 

i-5^/.50f ^ (203) 

i 

So far so good. What we need now is a local approximation for the two point function of An^^., 
the noise, which should also include the effects of the decay rate of the particles. The expression 
is given by [271] 



(An,.(t,k)An,,(t',k')) ^ (2vr)3j(3)(k + ^^ Y"^^^/,^^ 0{t - ufj^e{t' - t,) 



a^{t') ' 



(204) 



This kernel sources 5(f) only through the integral in cosmic time of the Green's function. We can 
then approximate this equation by (see Eq. 3.20 in [27,] ) 



(An,(t,k)An,(i',k')> ^ (2vr)35(=^) (k + k')^^^^ ^ .,(2vr)35(3) (k + k')^^, (205) 



with 



its 



2r^ 



(206) 



and N^iis ~ (f>/{HA). The computation of the power spectrum then follows as in sec. EJ 



Finally let us add a few words on our assumption that is a light field. If they were too light 
one might wonder whether they could contribute to the density fluctuations. Notice that since 



Notice that in generalizing the result derived for the background into the one for the perturbations, we are assuming 
that effects due to the extrinsic curvature of the surfaces of constant (f} are suppressed by powers of <C 1 i.e. 
7 ^ g\(f>\/H , at the time of freeze-out with k ~ {g\(f>\)^^'^ . 
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we require H <C F^. < M^. (r^^.), we do not necessarily need to be as light as Hubble, in 
fact rriy, < M■^^{T^^) /2 would be just fine. In this case they quickly redshift away. (Let us stress 
that even if they were effectively massless they could still be irrelevant for the late time curvature 
perturbations.) 



8.1.2. Non-Gaussianities 



To compute non-Gaussianities we follow [27\ and replace (p hy (p + 5(j)i + 5(j)2 in Eq- (|197p . and 
expand to linear and quadratic order in 54>2 and 5(pi respectively. To this purpose we need to 
compute particle creation in a time dependent background to second order in the perturbation, 
namely ^^^^n^. In Eq. (jl99p we neglected the space variation of the field in the computation of Sn-,^, 
but at the non-linear level those gradient terms will provide us with the largest non-Gaussianities. 
(Recall the term proportional to 7(9j7r)^ in Eq. (jl36p is the one that induces the largest effects 
since A;^ ~ ^/^H.) 

Given that the Oj = Xi are scalars, as we argued in sec. 17.31 we are then able to extend the 
result in Eq. (I194p to a spatially varying field, namely (j) ~^ n^d^(j) with ~ df^4>, hence 



5/2 

. (207) 



r^|A|(2^)3 r^|A|(2vr)3 " 

(This is analogous to the covariant version of Schwinger pair production in an electromagnetic field 
with a purely electric background E ^ <f), which turns out to be a scalar function of = F^^F^" 
and B = F^.F;,.) 

One might worry about the validity of this procedure once we add inhomogeneities. In our 
example the particles are created in the state given by Eq. (jl88p . where (momentum) gradients 
of the unperturbed distribution are suppressed by factors of 1/ {g\4>\Y^'^ = In fact, k is the 

scale that controls the validity of the local approximation. If we ignore these higher derivative 
(extrinsic curvature) terms, which are suppressed by /c*/k, we can then choose coordinates 
such that i = t + IT. Hence for equal i surfaces we get (^df:(p{t)^ given by Eq. (jl89p . Then 
we just replace — )• n^d^ for any coordinate system, with ~ g^^dyt. To incorporate all 
type of corrections we also need to include terms involving the extrinsic curvature, 5Kf^^, and so on. 



From the expression in (I207p we already obtain the terms (with vr = 

if 2 - cldfn2 + i^H + 7)7^2 + 7 (a^l - lidiTri)^ + ... = Noise, (208) 



with a some numerical coefficients (which may vanish). In particular we get 7(5j7ri)^ at the 
non- linear level as advertised^^ , leading to |/nl| — ^ 1 for 7 ^ ff, similarly to the calculation 
in sec. 17.11 



There are yet other sources of non-linearities which can be induced and we have not incorporated 
(represented by the ellipses). Since fi ~ cp'^, one may worry about terms like fiirdOi, which do 
not respect the shift symmetry. (Moreover, they could induce non-negligible non-Gaussianities in 
the squeezed limit). However it turns out these terms actually cancel out, and the theory has an 



Notice that for the model at hand we also need to include the perturbations in (</>). This will slightly modify 
the form of the function of \J —dc^P- in Eq. (|207|l . and ultimately the coefficient a in Eq. ((208]). 
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approximate shift symmetry, as we found already at linear order. To show this let us return to the 
interacting part of the effective action in Eq. (jlBip 

and perform the shift (p ^ (p + c. Assuming V{(j)) obeys the slow roll conditions, the action will 
remain invariant provided we shift at the same time (pj — >• = — c. The new action then has 
the same form as the original one, with — )• (pj. Since, as we saw previously, the strength of the 
trapping mechanism depends on the velocity as it passes through the sweet spot, the physics stays 
essentially unchanged as long as we sum over a large number of periods. Hence the EOM will be 
(approximately) shift invariant. The crucial point is that the change (p ^ i;^>-|- c is compensated by 
a concurrent shift in ti, i.e. ti ^ ti + 6ti with 6ti ~ c/(p, which is ultimately summed over a large 



number of periods [27[. Although this is not entirely obvious, the would-be breaking terms already 
canceled out at linear order. 

Something similar occurs in the model of [27| before making the local approximation. In that 



case the contribution to the EOM from particle production reads |27l | 

rh^6^ + y* dt'm^ (^^H{t') - m5cp{t')^ (210) 



which appears to have a 'mass' term and violate the shift symmetry. However, that is not the case 
since the change in that term cancels against the shift in the second term in the integral, using 
a{t) ~ e^*. In the local approximations these pieces do not even show up, since they explicitly 
cancel each other out. To include these effects to all orders we need to be careful with the sum 
over fiOi. At the end of the day the same cancellation occurs at second order and so on, such that 
we restore the shift symmetry as we argued before. 

Finally let us stress that our results here apply for the case of a local response. In principle. 



as in the trapped inflation scenario of 27|, large non-Gaussianities may also be produced via the 
non-locality in time of the response. 



8.2. Warm inflation 



Another example where we have dissipation/fluctuation is warm inflation [3414361] . In this case 
O represents ADOF in thermal equilibrium at a given temperature T mutually interacting with 
the inflaton (p. Even though Tq^ may be large enough to modify the dynamics of (p^ it remains 
sub-dominant with respect to the potential energy, V ~ MpH'^, and one could hope inflation might 
not be drastically perturbed as long as T <C V^^'^. The interesting aspect of this approach is the 
possibility of having T ^ H, which implies that thermal fluctuations dominate over vacuum effects. 
In warm inflation the evolution of the inflaton is governed by an equation of the type ^] 

4>+{3H + j)^ + V'{(P) = K^, (211) 

where ^ represents the (Gaussian) noise with (^) = and K ~ V^fT. Provided 7 does not depend 
significantly on the temperature, the computation of the power spectrum follows the same steps as 
in sec. [H The cas e dry ^ is more elaborate since one has to include perturbations in T already 



at first order (see 46|, |47| for more details). 

The dependence of 7 on the temperature varies with the different realizations of warm inflation. 
The most promising example is given by the so called two-stage decay model [36| . Similarly to the 
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trapped inflationary case (although due to different mechanisms) the rolhng inflaton produces some 
(ultimately heavy) particles (m^ S> T) through a g^^^x^ coupling, which subsequently decay into 
lighter degrees of freedom, t/j's. As we discussed before, this 'catalyzing' mechanism may allow us 



to approximate the evolution equations by local dynamics. However, contrary to the model in 271], 
in warm inflation the decaying product is assumed to have thermalized. The friction coefficient 
can be then computed and scales like [s^ 

jiT)^g'h^(^Y^, (212) 

where h and m = gcp enter in the coupling between x ^^'^ y^i G-g- ^xv ~ hmx^'^-^^ Clearly 
this effect is rather inefficient, since rriy. ^ T. One possibility is to increase the number of x 
particles such that more than one field is excited, or equivalently the number of decaying channels. 
Unfortunately, if we denote by M this 'enhancement factor', the assumption that warm inflation 
occurs for sufficient e-foldings (N^ > 50 - 60) with 7(r) > H requires M > 10^^,^.^^ 



Non-Gaussianities in the warm inflationary scenario have been also computed in the literature. 
In most cases the contributions are small, i.e. slow roll suppressed [HS] (since they stem off the non- 
linearities induced by the potential). Following the reasoning of our previous section, to include 
other type of non-linearities in warm inflation one needs to generalize the 7(/) coupling to the case 
where — )• -|- 6(f), which we would now write as "fu^d^cp with u'^ some four vector. In principle 
there are different possibilities for depending on the model. For the two-stage case we can have 
either (as before) or u^-, the four velocity of the x particles. The latter indeed would enter 
through the computation of the 7 factor in Eq. ()212p . However, in the limit T <^ m^, the x 
particles are non-relativistic (namely they are created essentially at rest) hence ~ n'*.'^^ This 
suggests a coupling of the form 'j{T)n^dfj_(j) with ^(T) given by Eq. ()212p . which following our 
previous arguments, would lead to non-linear effects of order |/nl| — l{T)/H. 



9. CONCLUSIONS 



In this paper we gene ralized the EFT for single field inflation [13h23|| to include dissipative 
effects. Following 13, 31 1 we introduced ADOF described by composite operators in the effective 
action compatible with the symmetries of the long distance physics. Our general assumptions were: 
i) the existence of a preferred clock, ii) negligible (or vanishing) contribution from the ADOF to 
the curvature perturbations at late time, iii) the time scale for dissipation and fluctuation induced 
by the ADOF is much smaller than a Hubble time with negligible memory effects, and iv) the 



In |48l ] it was argued that assuming warm inflation occurs and the potential obeys a series of slow roll conditions [36| , 
the thermal hypothesis becomes plausible and prad — {j/H)^'^ has a stable equilibrium provided j{T) ~ T^, as in 
Eq. ((2T2)) . 

Even though one may be able to relax the condition on the curvature of the potential (the so called 77- problem), 
some sort of tuning reappears in the form of very peculiar conditions on the matter content of the theory. Moreover, 
another important challenge for warm inflation model building is to keep under control the radiative and thermal 
corrections to the effective potential which could potentially ruin the slow roll conditions. (More so if one is going 
to assume 7 becomes sufficiently larg e [Hit , in view of the previous requirement.) Supersymmetry may be invoked 
to tame the radiative corrections [52[ ■ However, Supersymmetry is broken at finite temperature (and for non-zero 
vacuum energy), therefore some extra tuning may be required. A detailed account of the tuning of thermal inflation 
models lies beyond the scope our present paper. 

This is a consequence of the validity of a derivative expansion for heavy m^'s, which also supports the local 
approximations. The replacement u'^ — > n'^ agrees with the expressions in [i?]] . 
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validity of an expansion in spatial derivatives. The last two allow us to use a local approximation 
for the dynamics of the perturbations. Under these conditions then we showed that in the strong 
dissipative regime the two point function is dominated by the noise, and that it can be significantly 
larger than the standard result for slow roll single field inflation. In particular 



2c* (ctNc 



,2 ' 



(213) 



which departs considerably from the Bunch-Davies result. The reason is twofold. First of all the 
size of the perturbations for the canonically normalized vr field (tTc = ^/IV^tt) is larger than in the 
Bunch-Davies vacuum, but moreover because the normalization scale Nc can be smaller than the 
value it takes in single field infiation without ADOF for fixed H and Cg, namely c^N^. < 2Mp\I 
This is even more transparent if we assume the FD theorem applies, or formally define To as 
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in which case 



ToH'i 



2ct 



CshTpHj 



(214) 



with 

2M^\H\cs >Af,we get^^ 



c^Nc, similarly to what happens in warm inflation 34l436l|. Then, introducing A^ = 



iCCh 
(CO 



BD 



Ar 



»1, 



(215) 



already for To > H when j ^ H. The factor of Cgk^To in the numerator can be understood as 
follows. The energy density for vr is given by {diTc)^ ~ ^IttI (where d represents the derivative 
with respect to = x"^ /cs), which goes like k^To (a.k.a. Rayleigh-Jeans law). Hence vr^ ~ k^To = 
Cgki^To, which we compare with the quantum noise in the Bunch-Davies vaccum given by H^. 

As a result it is no longer the case that the power spectrum provides us the values of H^, 
and as in the standard scenario, but rather with a set of new parameters: ^,vo,Nc (or 
To = T^o / {^d))- Therefore, the tensor/scalar ratio would no longer give us a clean measurement 
of the energy scale of inflation. (Also because the ADOF might significantly contribute to 
gravitational wave production [s^].) 



We were also able to identify some specific signatures for non-Gaussianities. In particular, we 
showed that the level of non-Gaussianities can potentially be much larger than in the case of single 
field infiation without ADOF by a factor of ^/H ^ 1. Our main observation was the following: 
in order to dissipate we require (/>'s in the EOM, however, we need not only perturbe (j) but also 
'the dot', since the equal time surfaces are fiuctuating. Hence, the non-linear realization of time 



Notice this implies that the cutofT scale for higher dimensional operators in the EFT, A.u , (such as those obtained 
from (1 + (7°")" type of terms) may be lower than in the case of single field inflation without ADOF (for a given 
value of H,Cs), where Au ~ Mp\H\Cs for Cs <C 1 [HI. However, keep in mind the our EFT description breaks 
down when E > Fo, hence the 'lowering' of Au is only meaningful in cases where Au ^ Fq. 

One might be tempted to identify c^Nc with the symmetry breaking scale Af [2l|], but this would not be correct. 
The reason is that in order to define the scale at which time translations are broken one has to properly account 
for the contributions to T'^'^ stemming from all degrees of freedom, including the Cs. Following the procedure 
outlined in [23|, one can easily show that remains at 2Mp\H\cs- (This is actually not that surprising, since the 
normalization of the Goldstone boson it did not change, as show in Eq. (ff8|l.) 
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diffeomorphisms requires the following combination 



F{dt^) ^ ^ F((/)V^), (216) 

which not only induces 77T dissipation, but also (among others) a non-linear term: —^(diir)'^. 
(Notice the relative sign is dictated by the non- linear realization of the symmetry.) Since horizon 
crossing happens at Cgki, ~ VlH, this type of non-linear coupling leads to 

^^&^~/nlC^I/nl|=^^. (217) 

The shape is plotted in Fig. [2l and peaks at the equilateral configuration. (This is not surprising 
given the fact that the non-linearities involve derivatives oi C,.) However, there is also a significant 
contribution at folded triangles xi = 1,X2 — 2:3 ~ 1/2. Other non-linear terms may depend 
explicitly on the noise, such as 50st^, and are plotted in Fig. [31 Despite the fact that it scales 
with a single power of tt one can show that its contribution in the limit ki ^ ks is suppressed by 



(ki/ks) with respect to the local shape, in agreement with the results in 451]. We will analyze 



the squeezed limit and consistency conditions in the presence of dissipation in future work. 



In this paper we also studied specific realizations of the type of operators introduced in the 
EFT and the matching procedure. In particular we analyzed a local version of trapped infiation 
where the produced particles decay after they are created, which leads to (approximately) localized 
response functions. We showed how the term 7(9j7r)^ gets generated, with the subsequent j/{c^H) 
imprint on |/nl|- Crucial aspects of the model include: i) The ADOF responsible for dissipation 
do not contribute to the density perturbations at late time, ii) The emergence of a shift symmetry 
at the level of the perturbations, and in) The response functions were predominately sensitive to 
the preferred clock (p, whose fluctuations uniquely control the end of inflation, via ~ 5^(t + vr). 

Intuitively, the necessary gradients of vr appear as a result of the fluctuations of the clock, the 
field (j) itself, which sets the equal time surfaces where the unperturbed computation is assumed 
to hold to a good approximation. The derivative expansion remains valid as long as the typical 
length scale for the variation of the extrinsic curvature (of equal time surfaces) is larger than the 
typical wavelength of the produced particles, namely 1/k <C l/k^,. (Our conclusions also apply 
to the two-stage model of warm inflation, provided one succeeds in producing sufficient e-foldings 
while having 7 ^ in a consistent fashion.) 

One might wonder about the possibility of having a 'second clock' controlling the response 
functions for the ADOF. As long as we are only concerned about effects on the dynamics of the 
one clock driving infiation (assuming this second clock produces negligible direct contributions to 
C), one can incorporate its presence in the C-system by replacing n'^dn — )• UQda- We will study 
this in more detail in future work. 



In a nutshell, departing from the vanilla single field scenario opens new possibilities which may 
well be realized in nature. Once again, the EFT machinery is a wonderful tool to reduce the plethora 
of conceivable realizations to a theory of low energy degrees of freedom coupled to a set of composite 
operators whose correlation function encode all the information about the dissipation/fiuctuation 
properties of each specific model. In our case we reduced the number of additional parameters to 
three: ^^vqi^c- (Also Fq and Mq, controlling the validity of local approximations.) By taking 
the ratio between (^-correlation functions, such as the two and three point functions, we manage 
to cancel out most of our ignorance on the underlying dissipative mechanism, thanks to the link 
between different n-point functions induced by the symmetries. In this fashion we were able to 
show (assuming the noise is Gaussian) that the bispectrum peaks at equilateral configurations and 
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moreover ~ ^3/7 ^ ^^^^ class of models. 

A detection of the generic type of signatures we discussed in this paper, as a result of incorpo- 
rating dissipative effects during inflation, would increase our understanding of the dynamics of the 
early universe and also lead us towards a more precise description of the inflationary epoch. 
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Appendix A: Examples with Ohmic behavior 



Let us consider the well-known case in which the system (described by vr) is linearly coupled to 
a "bath" of harmonic oscillators, such that 0,133 



Aot = >C^(7r,7r) +£int(g,7r) + Cq{q,q), 



(Al) 



with 



2 2 
Then the EOM become 



(A2) 

(A3) 
(A4) 

Plugging the solution to Eq. (IA4p back into Eq. \A2>\ we obtain an equation of the so called 
Langevin form: 



vr 



+ WqTT + ^ Caqa = 



vr 



where 



+ UjItI + J dt'^{t-t')TT{t') = J{t), 



J{t) = -Y,Caqr'it), 



(A5) 



(A6) 



with (t) representing the classical trajectories of the oscillators in the absence of vr and 



(A7) 
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In the above expression we introduced 

G°.(t-0 = '>(i-tO ""'"°""''» , (A8) 

the retarded Green's function for an harmonic oscihator of frequency cJq,. If we now assume 
random initial conditions, taken independently for each oscillator in the environment, we can 
consider J{t) as a stochastic variable. Then, as tt is being affected by the 'noise' induced by J(t), 
its dynamics acquires a stochastic character. In the jargon of particle physics, we have integrated 
out the bath of oscillators and obtained an effective EOM for the vr field, with a generalized 
friction term plus noise. (However, as we emphasized in this paper, this EOM does not derive 
from a Lagrangian of the form £(7r, vr).) 

As it stands the EOM for vr is non-local (unless very particular properties for the Cq, coefficients 
are assumed 2^, 37]). Nevertheless, as we argue below, there are situations where we get 



Im7(cL!) ~ 70;, (A9) 
with 7 a constant, such that the EOM takes the desired form 

vf + 77T + (Jovr = J (AlO) 
with (J) = (and after we redefine ujq to include a renormalization piece). 

A Lagrangian description of the sort of Eq. ()Aip arises in the so called Rubin's model, where 
a heavy mass M is coupled to a half-infinite chain of harmonic oscillators of mass m and spring 
constant mufj^/A, after diagonalization {ujr is the highest attainable frequency) In such case 



one has [37| 



Im7R ~ u; ^^^Y^ 1 - '^V'^l ^i^R - (AH) 
which is of the type in Eq. ([2]) for w <C ojr with 7 ~ {m / M)uj r.'^^ 

The reason non-local effects make appearance in the EOM relies on the ability of the envi- 
ronment to back react in our system. Hence in order to have a local approximation, we should 
consider a situation where the energy of the system is damped into the environment and the latter 
has a negligible back reaction effect (or it becomes important on a time scale much longer than 
the length of the experiment). 

As an example, and in the spirit of Rubin's model, let us take the case of a ring of mass M 
attached to a half-infinite rope. If we denote by y{x, t) the height of the rope as a function of x > 
and time, and place the ring at 7r(t) = y{0,t), the EOM read 

(A12) 
(A13) 

x=Q 





d^y 










Mix 




dy_ 






dx 



*° Notice the asymptotic beliavior is given by yR{t) ^ sinujRt for t ^ <^ii' , wliicli is milder than in Eq. ((Tj. 
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where F is the tension per unit of length on the rope (and we work in units where the speed of 
propagation is taken to be one). From the wave equation we know 




(A14) 



and choosing boundary conditions such that only outgoing waves are allowed (i.e. y{x, t) = f{x—t)) 
we obtain 

if + 771- = 0, (A15) 

with 7 = F/M. We can now simply add a spring of frequency ujq attached to the ring to return 
to Eq. (fHOjl . in this case with J = 0. (Notice that the rope represents the continuum limit of 
Rubin's model, where we take m — )• and ujr — )• 00 while keeping mujR finite, so that Eq. (jAlip 
leads to a constant for all times.) 

In a more realistic setting we may imaging fixing the rope (now of length L) at an end, such 
that waves will bounce back on a time scale of order ts ^ L, introducing non-local effects. In such 
scenario, and as long as we are interested in time scales t <^ ts, our local approximation remains 
valid. 



Appendix B: The optical theorem 



Even though in this paper we deal with dissipation, our results are still consistent with unitarity. 
To make the connection more transparent let us consider the forward scattering amplitude for the 
TT particles, which we depict in Fig. [5] as a 'self-energy' diagram, with vr propagators represented 
by the wavy lines. 





FIG. 5: The wavy lines correspond to tt and the dark interactions are insertions of On couplings. The line 
connecting the dots represents Feynman's time order product of Eq. (jBll) . 



As it turns out, this amplitude is proportional to the time order product 

{T{60ix)60iy))), 



(Bl) 



also known as Feynman propagator, iG-p{x — y), which appears after we integrate out 60 in the 
usual path integral formalism with Ci^t = SOtt. Unitarity, in the form of the optical theorem, tells 
us that the imaginary part of this amplitude must be related to the total power loss 30|,l3]|- Hence, 
using 



ImAn^n — ImGF(w) 



(B2) 
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and (for oj > 0) 

ImGF(a;) =ImG^t(w), (B3) 

multiplying by a factor of oj (to go from rates to energies) we obtain 

dE 

— ~ wim A^^^ ~ wlmG^^t (a;) . (B4) 
If we require this expression to match the condition 

dE n 

that follows from a dissipative term of the form 77r, we conclude (for uj > 0) 

ImG?,t(a;)^7a; (B6) 
as expected. In general the expression in Eq. (jBSP may be more elaborate, but the procedure 



generalize to any function of ui |30|, ISlI . 



The reader may be puzzled about the appearance of the Feynman propagator rather than the 
retarded propagator. However Feynman boundary conditions can be used to relate the imaginary 
part of self-energy diagrams with the total radiated power (or energy loss), similarly to what is done 



in the EFT for gravitational radiation of [SHa, lUl • The crucial point is that the boundary conditions 
are so chosen to ensure 'in' and 'out' vacuum states for the (5Cs (e.g. no external gravitational 
radiation), but at the same time producing an imaginary part which precisely account for the total 
radiated power. Another way to see this is to notice that the total rate induced by the coupling 
SOiT is proportional to (ignoring factors of k for the sake of argument) 

J2{Tr,0\7r6O{t)\0,N){Tr,0\7r6O{0)\0,Ny = Y,{0\^C>{0)\N){N\6O{t)\0) = {0\6O{0)6O{t)\0) , 

N N 

(B7) 

which is nothing but the sum over the square of the emission amplitudes Att^n, for all N possible 
(intermediate) states. Then, using the relation (valid for a; > 0) 

J dte''^\0\6O{0)dOit)\0) =2Imi J dt e'''\0\T{5Oi0)5Oit))\0) , (B8) 

we reproduce our previous result (after multiplying by a factor of oo). This is nothing but the optical 
theorem at work. What turns out to be a bit more subtle is how to obtain the correct retarded 
boundary conditions of Eq. (I14|) from the path integral approach. However, this is possible in 
the so called 'in-in' formalism [2^. 37 . [38| . (See (4ol ] for a discussion in the case of gravitational 
radiation reaction.) 



Appendix C: Retarded Green's function 

In this appendix we provide some basic features of retarded Green's functions, and in particular 
we discuss an example where ImG^^ ~ w as we used throughout the paper. 

Let us start with some axiomatic properties for the retarded Green's functions and let us work 
at zero spatial momentum. First of all, from causality we know G^t(a;) is an analytic function of 
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u for Imuj > 0. Moreover, the imaginary part is odd in uj: 

ImGg,{-io) = -ImG^d^). (CI) 
Also, the real and imaginary parts are related via Kramers-Kronig relations, 

J-oo Uj' -U 

. _P ^^5£«&M, ,03) 



where P stands for the principal value. This implies 

dJ — - < oo, (C4) 

-oo 

from which we obtain 

ImGf^t(w ^ 0) ^ 0. (C5) 

Notice that Eq. (jCSp precludes the existence of a pole at a; = 0, however, it does not rule out the 
behavior of Eq. ()20p away from the origin, as we show in appendix [Hi 

Let us now study the example of electric conductivity. As it is well known the relationship 
between the current and electric field in a conductor is given by 

je = -iuja{uj)A{uj), (C6) 

with E(aj) = —iu}A{uj). Then using linear response theory one can show (Kubo formula) 

Rea(a;) = ^"^^-^("\ (C7) 

UJ 

where G^^^ is the retarded Green's function for the electric current jg. In general this Green's 
function can be parameterized as 55|] 

with Mj{uj) a 'memory function', and aj a constant. Depending on the system, for some range of 
frequencies one can approximate Mj{uj) ~ i/tj + 0{ujTj) with tj the 'memory time', such that 

GIM:.--^ IrnGU^)^.^^, (C9) 

with ojd — l/'^j) = ctjUJD- For ^ we obtain the behavior as in Drude's model in Eq. ([7|). 

Notice that we also have Re Gl^^ ~ ctujd (which is obviously consistent with the dispersion 
relations of Eq. ()C2p ). This means, in principle, that we get a large correction in the real part 
of the Green's function. In cases where our O operators have a Green's function of this form, 
this could potentially lead to a large mass for vr. In our cases of interest we assumed there is a 
mechanism that forbids a large contribution such that itItj remains 0(e), as in cases where there is 
an approximate shift symmetry. 
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Appendix D: Decoupling limit 

In single field inflation one can show that 6N and 6Ni are suppressed in the slow roll 
approximation, so that one may work with the theory of vr up to 0(e) effects, for models not too 



dose to de Sitter [iSy. Here we argue that decoupling still occurs even after we include ADOF, 



and make some general comments about the structure of the mixing terms. In this section we 
work in Mp = 1 units unless otherwise noted. 

As we know, the full T^'^ for the background takes the perfect fluid form 

Tfj,u = {ptot + Ptot)Ufj.Uu + (jtiuPtot (Dl) 

(which follows from the isotropic and homogenous conditions). We will work in the Newtonian 
gauge where 

goo = -1 - 2$, go^ = 0, gtj = g^j{l - 2^). (D2) 

If we concentrate on first order scalar perturbations we have 

5Tij = -2^!gijptot + gij5p + 5Tij (D3) 
5Too = 2ptot^ + 5p (D4) 
^Tio = -{ptot + Ptot)dM, (D5) 

where we included a dissipative term Jr,;,-, induced by the ADOF. (Keep in mind that the ADOF 
also enter in 5p, 5p.) The EOM become |56| 

dj {6p + 2do{eH^5u) + 6eH^6u + 2eH'^^) + didr'j = (D6) 



1 



5p + 3H{5p + 5p) + Hdr'i + 2eH'^ { - 3^ ) =0 (D7) 



-Sp + 3€H^6u = (D8) 



a2 2 

^ (d,d, - ls,,V^^ {^-^) = (<)r; - ^5,, Jr/) (D9) 

where we used ptot + Ptot = —2H = 2H^e, (e = —H/H^). Let us assume for the moment that 
5p, 5p, 5u and Srij do not depend on the metric perturbations. Then, from these equations one 
can immediately note that the metric perturbations decouple in the e — )• limit. Moreover, the 
constraint equations (last two) give us the value of ^ and ^ in terms of 5p, 5u and Srij. At the 
end of the day we end up with equations where the mixing with gravity is suppressed by e, and 
moreover, is independent of the relationship between 5p and 5p. 

To analyze the relative importance of the contributions of these mixing terms, we can combine 
Eqs (lDi-[D8l) to obtain (A: 7^ 0): 

6h + ^H{5pk + 5pk) - ^ [dpk + \5t^ + ^ {H{5Ti)u) - 2eH^6pk - SeHHpk 
+2eH'^5Tk + WeH^—Su - 3—^ {6pk + 2H6pk + 4.H'^6pk) + ©(e^) = 0. (DIO) 

Qj K 
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In this expression we replaced di{6T'^)k — >• 2/3 dj5Tk, with St^ a scalar mode which follows from 

the decomposition ((5rj)J^ = {k'^y — l6ij)6Tk + ■ ■ ■ (P = 1), and we used Eqs. (ID81 IDOP to solve 
for <i>, 

$ = ^Q5pfc + 5rfe^ +0(e). (Dll) 

Therefore the presence of gravity has two net effects. First of all, there is a "mass" term 
proportional to \^H; and secondly there are mixing factors, of order eH'^/k'^ at horizon crossing 
(recall k^, = k/a{t-f,)). Hence, as long as cj* ^ e^^^H and k^, ^ e^^'^H, we can ignore the mixing 
with gravity. For us, since k-^ ~ ^^B. jcg > H, this suppression is larger than the usual case upon 
noticing 

mixing ^ < 1. (D12) 

So far we have assumed that terms involving the metric perturbations in 5p, 5p, 5u and 5Tij 
are negligible. To clarify this point let us take the example of single field inflation without ADOF 
[l^ . To linear order in the perturbations, the stress tensor obtained from the first contributions 
to the action defined in Eq. ([65]) . i.e.^^ 

\ j + M|(l + g00)2}, (D13) 

can be written as in Eq. (jPSp . using 

5p = ^7r + (p + p + 4M|)(7r-$), (D14a) 
5^ = ^vr + (p + p)(7T - $), (D14b) 
5u = -TT. (D14c) 

For these variables to be approximately independent of $ we need $ ^ tt. Solving for ^ using Eq. 
()Dlip without the ADOF we obtain (in the slowly varying approximation and restoring Mp) 

Therefore the assumption $ <C vr is self-consistent provided 

(, + , + 4M|)«l k.»^I±l±^^. (D16) 



Using p + p ^ 2eMpH'^ we recover an expression similar to ()D12p (but without the factor c^^), 

£ H 

mixing ~ < 1. (D17) 

7* 

Note that the above estimations do not apply in the limit e — >• 0, since the higher derivative 
contributions can no longer be ignored. To analyze this case we add the term proportional to M| 



Note that the term M|(1 + g°°f aUows for / 1 (see Eq. (fff)) ') 
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in Eq. 



I.e. 



Setting p + p ^ and — >■ we get 



ru2 



1.2 



;7r - 3{H<^ + 



(D18) 

(D19a) 
(D19b) 



(The second line is valid provided <C vr and + ^' <C ^tt.) Using again Eq. (|Dlip (without 
ADOF) we have 



A;2 



Sp 



4M|— vr 



(D20) 



Therefore we see the assumption that we can neglect the dependence of Sp on the metric pertur- 
bations requires, not only M2 ^ Mp, but also [13] 



fc2 



h > M^/Mp 



(D21) 



In addition, this explains the apparent puzzle in the mixing expression of Eq. ()D12p . that appears 
to exactly vanish in the de Sitter limit. (This we know is not the case, for example in the Ghost 
Condensate. See 13(] for more details.) Notice that the condition in ()D2ip does not make any 
reference to Cg, which is zero in the de Sitter limit. In general, for small (e, Cs) the mixing with 
gravity will depend on which one is larger between p + p and M2 , as in (ID16|) . and therefore in the 
near de Sitter limit the mixing may become enhanced with respect to e. 

The inclusion of ADOF does not modify the previous analysis considerably. Similar arguments 
can be made for vector and tensor modes. 



Appendix E: Comment on OSK^ 

Let us briefly analyze the importance of the terms with the extrinsic curvature. In particular 
we compare the operator s{t)^^SKli, given in Eq. ([89]) . with f{t)Oi of Eq. (f87|) . For the sake 

of argument we take Oi ~ O. By this we mean that for these operators all response and noise 
coefficients are of the same order. Then to linear order in vr, from f{t)Oi we get firOi, whereas 
from the one with the extrinsic curvature we have ^rfi-O^n. 

Taking the ratio we see that if A;^ satisfies 

A;2« MKf{t)/s{t), (El) 

we may be allowed to neglect operators with SK. As we have shown in this paper, the computation 
of the power spectrum and non-Gaussianities (due to the behavior of the retarded Green's function 
G^) is dominated by values for which, at freeze out, /c^ ~ yj / Cg ■ Therefore, to satisfy the 
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condition in Eq. (jEip we need 



<. Mk. 



(E2) 



Otherwise, the contributions from OSKji could be important. 

Appendix F: Dissipative Green's functions in an expanding universe 

In this appendix we discuss some properties of the Green's function relevant for the computation 
of non-Gaussianities in sec. [TJ in particular G'y{x,y) = y'^g^{x,y) {x = —kcgrj, y = —kcgT]'). In 
Fig. [6]we show some graphs of 6*^(0, y) for three different values of 7. (For a fixed x < (4+7///)^/^, 

the behavior of Gry(x,y) is similar to the one shown for x = 0.) We notice that peaks at 
y ~ (4 + 7/if)-^/^, and as 7 increases its amplitude decreases, whereas its support increases. 




FIG. 6: The Green's function G^{0,y) 
A + j/H^7^ (dashed) and 4 + -f/H 



y^g^{0,y) for three different values of 7: 4 + j/H — 5^ (black), 
10^ (dotted), the first peak of G-y is around y = 5, 7, and 10, 



respectively, supporting the fact that the location of the peak is at ^74 + j/H. 



For a particular value of 7, in Fig. [7| we plot the Green's function G-y{x,y) and its temporal 
derivative per Hubble time, dtG^{x,y) / H = —xdxG^{x,y), for x = y^4 + 7/// as a function of 
y > X, and for y = 1^4 + j/H as a function of x < y. To estimate the amplitude of integrals of 
it is useful to note that 



/" 

Jx 



dy Gy{x,y) 



-1. 



(Fl) 



Then, for large values of 7 a rough estimation for can be obtained by counting each power 
of as 7, adding a factor of 7"-^/-^ for each G^ (i.e. 7"^/^ for each g^) and a factor of 7-^/^ for 
each integral, with units made up by 'j/H. Note that having a time derivative per Hubble time 
of G^ cannot increase significantly (see Fig. [7|). On the other hand, the contribution of an 
additional spatial derivative per Hubble length will increase /^^ by approximately a factor of 7/c^ 



(since k'/{a'H') = y'/ci 



We keep the factor of 4 to recover the standard result in the case 7 — )■ 0. 
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FIG. 7: The Green's function Gj{x,y) (on the top) and its temporal derivative, per Hubble time, 

dtG-y(x,y)/ H = —xdxGj(x,y) (on the bottom); for x* = ^4 + 'y/H as a fimction of y > x* (on the 
left), and for y* = y/4~+j/H as a function oi x < y* (on the right). All plots correspond to x* = y* = 10. 



Appendix G: Mixing 



To study the effects of friction at leading order in tt, in this paper we concentrated on the 
study of a generic type of operator, i.e. Ovr, with O including in principle a series of contributions 
(after integration by parts). However, once we start adding higher order effects, we treated each of 
them separately. The reasons were both the fact that not all the terms generate the same type of 
non-linearities, and also obviously for simplicity. In general however more than one operator will 
be present at the same time. Here we make a few comments about the kind of effects that may 
appear as a result of including more than one type of ADOF concurrently. 



1. Scalars 



Let us start with scalar operators Oa, A = \,. . .N. In general, in addition to the contributions 

of the response and noise given by the self-correlation functions, we will also have the ones given 
by the mixed-correlations (e.g. {50^^50^)). To be more precise, we are considering interaction 
terms of the form 

Sint = - j d^xa^{t)OAF^, (Gl) 

where are taken as external forces that slightly disturb the dynamics of the ADOF associated 
to the O's, and where repeated index are summed over. The linear response of Oa due to the 
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application of the external forces can therefore be written as 

SO^ = - j d^x'a^{t')G%{x,x')F^{x'), (G2) 

and 

G%ix,x') = ie{t-t'){[60A{x),60B{x')]). (G3) 

Additionally, we have the fluctuations of the ADOF, which are characterized by the two-point 
function of stochastic sources 60j^{x): 

NSj,{x,x') = {50i{x)S0^s{^')). (G4) 

By concentrating in the case where a local approximation applies, as described in sec. I5.lt we 
can write 

60^{t, k) ~ k) + r^^(t) J-at(a3/2(t)F^(t, k)), (05) 

{50^{tM)S0UtM)) = {2Trf6^^Hki+k2)N^B{t,t'), (G6) 



with 



NUt,t')^ '^ABit)^-^^j^, (G7) 



where fiABit), ^Asit) and I'ABit) are approximately time-independent. 



00 



2. fit)Oi and O29' 

For the sake of completeness, let us study the effects due to the mixing between Oi and O2, 
under the assumptions of the local approximation described previously. Note that if we ignore that 
mixing, the interaction /lOi yields a term 771" in the EOM for vr (see sec. 17. 2p . while 02g^^ produces 
a higher derivative contribution if we assume that both operators have a Green's function satisfying 
ImG(a;) oc a; as w — )• (see sec. [73]). Hence it is reasonable to consider the case where the term 
/lOi dominates over f25g^^02 at linear order. This allows us to concentrate in a situation in which, 
at leading order in the corrections introduced by O2 , we can neglect the contributions coming from 
the interaction 6g^^02 alone, and its effects enter only through the mixing with Oi. For example, 
this would be the case if Oi = O2 and ^ /i- So in practice we can work at leading order in /2. 
Moreover, for simplicity, we concentrate here only in the contributions obtained using the linear 
response approximation, although as we have seen in sec. I7.5l the non- linear response contributions 
obtained from /lOi could produce the largest non-Gaussianities in this kind of models^^. 

There are two types of mixed response terms: one is given by how /i(t) (the force coupled to 
Oi) affects 5O2 and the other by how f^Sg^^ affects 6O1. Using the local approximation, these can 



''^ Note that from the analysis of sec. 17.21 if we have only the interaction /lOi the level of non-Gaussianities obtained 
within the linear response approximation are negligible. 
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be written as 



c-.^R ■ „ ( 9 5,-7r9,-7r 



/2ri2 ( 27f + 2nn - + ^^^-^ ] , (G8a) 



^2 



50^-- A21 (^A^ + Ay J + r2i (/i^ + fi^n) , (G8b) 

where /i^is = fiAB + 3/ 2HTab and we assumed the existence of an emergent shift symmetry, as 
described in the main text. 

To hnear order in /2 the equation for tti becomes 

^1 + {3H + 7)^1 + cl^TT, = (^h60f - 2/2^^^^^) , (G9) 

where = N, + 2f2h{T2i - Tia)), = - 2f2N-^ hcl{T2i - Tia), and 7 = 

^c~^ f A^rii + 2/2A(3iiTi2 - /ii2))- Note that c^Nc = c^Nc (recah we are working to hnear 



order in /2). To ease notation we omit tildes in what fohows. The particular solution to the above 
equation is 

fikcs n 

T^i{k,r]) = ^ / dr]'g^{kcs\v\,kCs\v'\)SOf{'k,ri') 



and the resulting power spectrum can be written as 



P^ = p/2=0(i + 6ll^) (Gil) 



r Q . 2 

where Pp~ is the one obtained for /2 = 0, which is given in Eq. ()132p where fo = J^ii/i • The 
source for vr reads 

J = Ar-iA(2/2OTi2-/2/ii2)^^-2iv-V2A(ri2-r2i)^^ 
+ 2/2iv-Vir2i^^ + /2iVc"'A(/ii2 - 6OT2i)^2 + 2/2iv,"Vi(ri2 - 2r2i)^^ 

_ A,of _ y/jkpl _ ,f/A^ ^ 2/.«-*(^. ,G12) 

Note that if Oi = O2 the second term vanishes. 

As we have learned from our computations in sec. [71 the contribution of the non-linear terms 
with spatial derivatives become more important for generating non-Gaussianities than those with 
only time derivatives. Let us then analyze for simplicity only the first and the last term, which 
we expect to be among the leading ones. Notice that the calculations for the first source term are 
exactly the same as the ones we have already performed in sec. 17.11 in that case for the first term 
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on the RHS of Eq. (|137p . In fact, it can be written as 



with 



1 diirdiTT 
Jk = -sil —, (G13 



fi Til /iFn' ^ ' 



where we have si <C 1 consistently with our hnear expansion in /2 after taking, for simphcity, 
^AB ~ HfiAB- From eq. (|144p we obtain 

which could be still large for j ^ H, even when Cg ^ 1 and within the validity of our approxima- 
tion, namely si ^ 1. 

Similarly, for the last part of the source term we obtain 

, 2 



/nl = ^ V / dy'y'^{g^{0,y')f 



X 



10,,2-«V^vr^r(^+4 
3 c 

+00 







3 di ^ 

dzz^g^ {y',z)g^{0,z), (G16) 



where S2 = ^^y^j^ ^ 1) and again for simplicity we take all mass scales to be of the same order. 
After a numerical calculation, in the regime 'j/H S> 1, we can well approximate /nl by 



/n1^-3^^. (G17) 



similar to our previous case. 



Vectors &; Tensors 



Finally, as we discussed in sees. 14.21 and 14.31 there is a family of vector and tensor interactions 
terms that can be added. These can be written as 

5^ = - 1 d'^xV^ {h(t)Ol + f2{t)Of + h{t)OT + ...) (G18) 

- 5^ [si{t)0\ + S2{t)Of + s,{t)OT + ...) + ••• 

For the sake of simplicity (and notation), let us here ignore the factors fi(t)^s and Sj(i)'s, whose 
time dependence is in general slow roll suppressed. 
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After introducing vr, these terms become 

O? = 0'^df,{i+TT) = O? + O^d^^TT (G19a) 
Of = Of + 7r)a,(t + ^) (G19b) 

0000 ^ 0M^'P5^(i + ^) (G19c) 

etc. Adding all them up we get 

O? + 0^° + of + . . . = O + O'^a^vr + O^^a^vra.vr + O^^^a^vra^vrapTT . . . , (G20) 

with O^ "'^ made up from some combination of our original tensor operators. We can rearrange 
the second line of Eq. (jGlSP in a similar fashion. 

As an example, let us consider O^ = O^j^g^^, which produces the first vr-dependent term in 
Eq. ()G20p . Notice that by simple inspection, at linear order in tt, we can transform the analysis 
for this operator into the one for a scalar operator we studied throughout the paper. Indeed, after 
integrating by parts we obtain 

j d^x./^O'^d^ir j d^x {d^^Of") ir = j d^x^On, (G21) 

with O = —^^fJ,{^/—gO'^), an effective scalar operator. By construction, this interaction obeys 
the shift symmetry. Moreover, assuming it satisfies the hypothesis of sec. 15.11 it may induce 
a large friction term for 'y ^ H, and yet no 0(7r^) terms in the action. Naively this might 
suggest an absence of connection between friction and non-Gaussianities in this case, however as 
we emphasized throughout the paper, the non-linear response will induce non-linear couplings that 
indeed we expect to be large, leading to /nl — l/c^H. 



Appendix H: OSg'^'^ model(s) with 771 dissipation 

In this appendix we discuss an example of a Og^^ type of coupling with 771" dissipation. We 
argued in sec. 11.21 this requires some peculiar analytic structure, in particular (see Eq. ([20]) ) 

luiGfM ~ (HI) 

within a range of frequencies fio < to < Tq to be determined momentarily. Here we introduce 
one possible model where we couple the inflaton to a scalar field with strong dissipative dynamics 
induced by a second (hidden) sector that does not couple directly to vr. We will explore this set 
up in more detail elsewhere. (We drop the tildes from now on, and due to the plethora of 7's we 
use 7,r to identify the dissipative coefficient for vr.) 

The basic idea is to use the EFT of inflation of [l^ but, in the spirit to the model in [21], let 
us couple vr to a scalar field a via a term a6g^^ so that we generate a pan coupling with p having 
units of mass [21]. For its dynamics we write 

J^a = \ {d^<Tf - + aJ^, (H2) 
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where stands for the interaction between a and a dissipative sector described collectively by 
which does not couple directly to vr.^^ In essence this is a two-steps model where vr dissipate 
via a mixing term. This can also be described in the framework of [i^, except that in addition we 
have a sector which couples to a and induces (strong) dissipation. (We could for example put the 
^'-system in thermal equilibrium at a given temperature T-q/.) 

If we now compute the retarded Green's function for a we have (see Fig. [8]) 

where is the self-energy contribution from the ([J^, Ji/,]) correlator due to the aJ^n coupling (see 
below). We assume now that the imaginary part of this self-energy can be approximated as 

ImF.,(w)~7^^^, (H4) 

which is a more standard behavior (similar to Drude's model in Eq. ([7|) and behaves linearly in 
(jo for oj < A^, see appendix [Cj Using analyticity of F^ in the upper half plane we also have (via 
Kramers-Kronig relations) 

ReFvi>(w) ~ 7*A* ( ^2^\2 ) - 7*A*, (H5) 

for u; < A^, . Therefore 

p^ImG^,,(^,k) J^^'" , . (H6) 




FIG. 8: The wavy lines correspond to tt and the dark interactions are insertions of a-k couplings. The larger 
dark circle represents the Green's function for the ([J^, J^]) correlator from to the crjip coupling. 

Our quest is to find a regime where we get 7,7 vr^ for the EOM, or in other words the scaling of 
Eq. (j20p in Eq. (jHOp . That is indeed the case when 

_ k2 - ^ 7^Avt] ^ < 7|a;2. (H7) 

Naively this seems like a difficult task given the constraint u ^ Aijr, however we can always set 



In principle ^ could talk to the inflaton via gravitational interactions, but these are subleading. 
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/i to cancel out the large contribution, for example by tuning^^ 

7xtAvt -/u^ ~ F^. (H8) 

Moreover we also assume u)^,,k^ <^ 7.J, (which we show is self-consistent), so that under our working 
hypotheses we get 

p^ImG,-,,M:^^i, (H9) 



for < w, |k| ^ 7iii(^ A,!,), hence 



7- ~ — , (HIO) 

7* 



which gives 77r ~ for p ~ 7^,. As we showed in this paper, with a 7,r7i" dissipation term we have 
/cj, ~ \J^T^H, therefore the conditions which led to Eq. (IH9|) are satisfied, provided 'y,^ ^ H. 

To compute the power spectrum and non-Gaussianities we also need the noise part. This follows 
also from the (indirect) coupling between vr and Jiji via a. This can be shown by explicitly writing 
the EOM (ignoring the Hubble expansion for simplicity) 

Ti-fc + ZcphTrfc = -p&k, (Hll) 
(^k + (fcph + /^^)<7fc = + pTTk- (H12) 

As we argued above, in the regime we are interested in we can write (treating as a particular 
operator 60) 

J* = 4 + /g,V = 4 + 7*A*.-7**. (HIS) 

where we used Eqs. ()H4p and (jHSP and added as before a stochastic source term J^. Therefore, 
choosing fi according to Eq. (]H7p we have 



7^,0-fc ~ pTTk - J|, (H14) 
and plugging it back into Eq. (|Hlip we finally obtain 

Vffc + kl^TTk + JnTTk = Jn, (H15) 

with = (/9/7<i') J^. Hence we are left with an equation as in ([2]) and the computations through- 
out the paper follow. In particular we get non-Gaussianities of order /nl — 1it/{cIH) as in sec. I7.1i 

Notice we can also construct a similar model for a Ofj^g'^^ type of operator replacing a by a 
gauge coupling gA^d^j,!:. Then Jij, would play as similar role as the usual electromagnetic current 
(see appendix [C|) . 



It is not surprising we require some fine tuning while dealing with (unprotected) scalar fields. One possibility is to 
think of (T as a pseudo-Goldstone boson, or introduce a supersymmetric version, with the scale of breaking near 
Hubble. (In that sense we resemble the EFT for supersymmetric multi-field inflation of see also [2^.) We 
leave this open for future work. 
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Appendix I: Consistency of local trapped inflation 

Here we summarize some basic constraints that guaranteed the validity of our approximations. 
For the background equations we need 

\'4>\ « 3if|0| « (II) 

so that we have a constant velocity solution 



(12) 



9 

In addition, we assume that the energy density is dominated by the potential energy, 

= ^ + y(^) + ^ M^^n^^ ^ V{c^). (13) 

i 

The generalized slow roll parameter e = —H/H"^ is given by 
Then, e ^ 1 implies 

< y ((/.), (I5a) 
5^M^,n^, «F((/<), (I5b) 

i 

which are the conditions required to make the approximation in Eq. ([I3|) . To estimate the LHS of 
Eq. ()I5bp we replace the sum by an integral as above and find 

Gathering all pieces together we collect Eqs. (fTMl) . (fWl) . (fTMl) . (ITT]) . (lT5a|l . (jl5bl) . as weh as the 
conditions A^ph <C k and /cphA ^ |0|, plus » ff. In addition we also have the number of 
e- foldings, where we get 

N. = [^d<p = ^, (17) 



using (j) = (pi ^ (pE, with (pj (</>£;), the inflaton field at the beginning (end) of infiation. 
Next we analyze these constraints for the paradigmatic example 

V{cP) = ^<A'. (18) 
We have five parameters in the model: m. A, F^, g and the initial value of the inflaton fleld. 
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1. X 10" 



0.000 0.002 0.004 0.006 0.008 0.010 



6. X 10" 



4. X 10" 



3. xlO" 



2. X 10" 




0.00 0.02 0.04 „ 0.06 0.08 0.10 



FIG. 9: The region of the free parameters (m and g) allowed for consistency (with " ^ " replaced by 
" < "). The horizontal axis represents g and the vertical axis m/Mp. Top: = lOOi? (case I); Bottom: 
F^ = lOOOTJ (case II). 



6/. To reduce the parameter space we impose A'^e = 60 and 



C 



H 



10" 



A;3 



(19) 



from where we obtain the condition (see Eq. (jl33p . using (t>^) 



(110) 
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Then, 

1623/16315/16^3/4 4/rJA|,^ 1179648 X 1034 ^7rll5l3/2^9/2^ 

^ |A|~ ^ £ (111) 

For simplicity, we analyze two cases: = lOOH (case I) and = lOOO-ff (case II). In Fig. 
we plot the 2-dimensional region for the parameters m and g for which ah conditions are satisfied 
(with " <C " replaced by " < "). The horizontal axis represents g and the vertical axis m/Mp-. 
case I corresponds to the graph on the left and case II to the right. As an example, taking 
Fx = lOOF, m = 10-^^Mp and g = lO'^, we obtain if ~ 1.3 x lO'^^Mp, j/H ~ 135, (f> ~ 3.3Mp, 
6tc ^ 1.98xl0-4i^-^ |A| ~ 1.29xlO-SMp and e ~ 1/32. Also, for F^ = lOOOF, m = 5x10-^^ Mp 
and g = 10-\ we get H ~ 3.34 x IQ-^^Mp, -f/H ~ 539, 4) ~ 1.63Mp, 5tc ^ 1.4 x 10"^i^"^ 
|A| ~ 2.27 X 10"'^ Mp and e ~ 1/32. In all cases j/H <C g\(p\/H'^, which guarantees the validity 
of the local approximation A;* ^ k. As a final check let us estimate the size of the curvature 
perturbations one obtains from the fiuctuations in the additional scalar degrees of freedom we 
added into the theory. Following [22] we can estimate this contribution by computing 

K^^Sx - M^An^, (112) 



where we take M^{6tc) — g\4>\5tc — y g\4)\. In the above expression 5g^ is not the curvature 
perturbation (^y., but nonetheless it gives us an idea of the size of the curvature it induces. For the 
sake of comparison, let us evaluate this expression at k/a ^ H, where we get 

On the other hand, from the analysis in sec. [6] we have (see Eq. (j206p ) 



(C^C0> ^ ^/Hh^'-^^ —■ (114) 

Hence, using (j)^ ~ e^V (e^ < e), then 



i^9x^9x) Tx^ 



X^^2 



4 « 1' (115) 



since F^ < , /c^ < and ^ 1. Therefore curvature perturbations become subleading. 
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